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SECTION  I 


INTRODUCTION 

The  transonic  wing  design  procedure  (see  Fig.  l)  described  in  the  previous 
volumes  of  this  report  utilizes  the  following  computer  programs:  Lockheed's 
version  of  the  NASA  Ames/Bai ley-Bal lhaus  3A  extended  small  disturbance  (ESD) 
code1*2  (ATWP) ,  Jameson's  FL022  full  potential  equation  (FPE)  code3,  Vander- 
plaat's  constrained  function  minimization  (CONMIN)4*5  code,  McNally's  2D 
boundary  layer  code  (TWODBL)6,  a  wing-pylon/nacelle  mutual  interaction  code7, 
and  viscous  linking  codes  (LINKUP,  LINKDN).  Both  of  the  inviscid  wing  analysis 
codes  have  been  coupled  with  CONMIN  to  provide  a  set  of  programs  that  can  be 
used  to  design  and  analyze  wings  that  are  in  some  sense  optimized  for  a  given 
set  of  design  conditions  and  constratins.  The  boundary  layer  program,  TWODBL, 
is  used  to  compute  displacement  thicknesses  that  are,  in  turn,  used  to  compute 
both  "fluid"  and  "hard"  wing  ordinates  during  the  design  process.  The  wing- 
pylon/nacelle  code  is  used  to  evaluate  the  effects  of  wing-pylon/nacelle 
interactions  on  a  design  wing  pressure  distribution.  The  purpose  of  this 
User's  Guide  is  to  present  the  information  required  to  use  these  programs  in 
the  manner  specified  by  the  design  procedure. 

The  following  sections  present  general  descriptions  of  the  programs  and 
their  input  requirements,  detailed  descriptions  of  input  variables  and  input 
deck  formats,  job  control  requirements  and  descriptions  of  program  output.  In 
each  section,  the  inputs  have  been  separated,  as  much  as  possible,  into  blocks 
or  modules  according  to  their  function.  Sample  input  decks  are  given  to 
illustrate  how  these  blocks  are  joined  together  to  form  an  input  deck.  Ex¬ 
amples  of  printed  output  from  the  programs  are  presented  as  well  as  descrip¬ 
tions  of  output  directed  to  mass  storage  for  post-run  processing.  Finally, 
functional  descriptions  of  each  subroutine  referenced  by  the  programs  and  a 
list  of  the  FORTRAN  logical  units  used  by  the  programs  is  given  as  an  aid 
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for  future  program  modifications  by  the  user. 
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SECTION  II 


I NV I  SC  ID  ANALYSIS  AND  DESIGN  PROGRAMS 
1.  GENERAL  DESCRIPTION  OF  THE  PROGRAMS 
a .  Analysis  Codes 

The  two  analysis  codes  used  in  the  transonic  wing  design  procedure 
described  in  the  previous  volumes  of  this  report  compute  inviscid,  i rrota- 
tional  transonic  potential  flow  about  30  wing  configurations.  The  small  dis¬ 
turbance  code  is  subject  to  the  following  restrictions: 

(1)  The  wing  thickness-to-chord  ratio  is  small  and  the  leading  edge 
is  not  too  blunt. 

(2)  The  angle  of  attack  is  not  too  large. 

(3)  Embedded  shock  waves  are  weak  (M  <  1.3)- 

(A)  The  boundary  layer  is  negligibly  thin. 

(5)  The  wing  lies  in  a  plane. 

(6)  Wing  sweep  is  not  excessively  large  (A  <  50°). 

The  full  potential  equation  code  is  restricted  only  by  conditions  2,  3  and  A. 

The  salient  features  of  both  analysis  codes  are  described  in  the  following  sections. 
The  Ames  Qa i ley-Ba 1 1 haus  ESD  Code  (ATWP) 

The  Ba i ley-Bal 1 haus  3A  code  solves  an  extended  small  disturbance  approxi¬ 
mation  of  the  exact  governing  equation  for  potential  flow.  The  partial 
derivatives  in  the  governing  equation  are  replaced  by  finite  difference 
approximations  that  yield  a  set  of  nonlinear  algebraic  equations  that  are 
solved  by  Successive  Line  Relaxation. 

The  Bai ley-Bal lhaus  code  uses  an  embedded  computational  mesh  system  to 
accelerate  convergence  of  the  relaxation  process.  This  mesh  system  consists 
of  a  wing  oriented  fine  inner  mesh  embedded  in  a  coarse  Cartesian  outer 
mesh.  A  coordinate  transformation  is  used  to  map  the  wing  planform  into  a 
rectangle  in  the  computational  plane.  Grid  stretching  relations  are  used  to 
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cluster  mesh  points  in  the  regions  around  the  wing  surface  where  large 
gradients  in  the  flow  variables  occur. 

The  solution  process  includes  alternating  sweeps  of  the  crude  and  fine 
mesh  systems.  As  a  first  step,  the  embedded  fine  grid  is  relaxed  with  the 
potential  function  4>  fixed  on  the  fine  grid  outer  boundary  and  Neumann 
boundary  conditions  imposed  on  the  wing  surface.  After  sweeping  the  fine 
mesh,  the  circulation  distribution  is  computed  using  fine  mesh  potentials. 
The  circulation  distribution  is  then  used  in  (1)  the  far-field  expression 
of  Klunker8,  and  (2)  in  the  solution  of  the  Laplace  equation  in  the  down¬ 
stream  Trefftz  plane,  to  compute  the  potential  on  the  crude  grid  outer 
boundary. 

Next,  the  crude  grid  is  relaxed  using  the  updated  potentials  on  the 
outer  boundaries  and  Neumann  boundary  conditions  on  the  airfoil  surface. 

With  the  crude  grid  relaxed,  the  potentials  on  the  fine  grid  outer  boundary 
are  updated  (by  interpolating  the  crude  grid  potentials)  and  the  process 
repeated  until  both  the  crude  and  fine  grid  solutions  converge. 

The  Jameson  FPE  Code  (FL022) 

Jameson's  FL022  code  solves  the  full  potential  equation  for  inviscid 
transonic  potential  flow  by  finite  difference  techniques.  Unlike  the  ESD 
code,  the  FPE  code  uses  a  single  computational  mesh  system.  Shearing  and 
square  root  transformations  are  used  to  map  the  wing  surface  into  a  plane 
in  the  computational  space.  The  wing  surface  is,  in  a  sense,  "unwrapped" 
about  an  arbitrary  singular  line  just  inside  the  leading  edge  of  the  profile 
at  each  span  station.  The  resulting  coordinate  system  is  highly  non- 
orthogonal.  The  program  uses  a  remeshing  scheme  to  speed  up  convergence. 

The  program  is  started  with  a  relatively  coarse  grid.  This  grid  is  relaxed 
a  specified  number  of  iterations.  The  number  of  grid  lines  is  then  doubled 
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and  the  solution  is  restarted  using  values  from  the  old  grid  system  inter¬ 
polated  onto  the  new  grid  system.  This  process  can  be  repeated  at  least 
two  times  to  achieve  the  final  desired  mesh  system. 

The  finite  difference  approximation  of  the  governing  equation  utilizes 
the  rotated  difference  scheme  introduced  by  Jameson  which  adds  a  directional 
bias  in  the  direction  parallel  to  the  velocity  vector  at  points  where  the 
local  flow  is  supersonic.  This  scheme  helps  to  alleviate  problems  induced 
by  the  change  in  equation  type  (i.e.  elliptic  to  hyperbolic)  that  occurs 
when  the  flow  becomes  supersonic, 
b.  Optimization  Code 

Vanderplaats 1  CONMIN  code  is  used  to  perform  the  optimization  of  wing 
sections  during  the  design  process.  CONMIN  is  a  FORTRAN  program,  in  sub¬ 
routine  form,  for  the  solution  of  linear  or  non-linear  constrained  optimi¬ 
zation  problems.  The  user  must  provide  a  main  program  and  subroutines  to 
define  the  optimization  problem.  Problem  definition  consists  of  specifi¬ 
cation  of  the  quantity  to  be  minimized*  (the  object  function),  variables 
which  can  be  adjusted  to  attain  the  minimum  (decision  variables),  and  con¬ 
straints  which  the  solution  must  satisfy.  The  basic  optimization  algorithm 
used  in  CONMIN  is  the  Method  of  Feasible  Directions.  While  the  program  is 
intended  primarily  for  efficient  solution  of  constrained  problems,  uncon¬ 
strained  function  minimization  may  also  be  solved;  the  conjugate  direction 
method  of  Fletcher  and  Reeves  is  used  for  this  purpose.  The  general 

*For  the  purposes  of  this  manual,  minimization  and  optimization  are  one  and 
the  same  since  CONMIN  establishes  an  optimum  by  minimizing  the  user  specified 
object  function. 
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minimization  problem  is  to  find  the  values  for  the  set  of  variables  V(l) 

(the  decision  variables)  which  minimize  some  object  function,  OBJ,  subject 
to  a  set  of  constraints  G ( J )  <0  and  a  set  of  upper  and  lower  bounds  VUB(l) 
and  VLB(l),  respectively,  placed  on  the  decision  variables. 

CONMIN  has  been  coupled  with  both  analysis  codes,  in  the  optimization 
mode,  the  analysis  code  is  called  to  generate  the  pressure  distribution  used 
to  compute  the  objective  function.  CONMIN  is  called  to  evaluate  the 
objective  function  and  set  the  values  of  the  decision  variables.  These 
decision  variables  are  used  with  a  set  of  shape  perturbation  functions  to 
alter  the  shape  of  a  desired  control  station.  The  analysis  code  is  then 
executed  with  the  new  wing  geometry  to  update  the  objective  function. 

This  cycle  is  repeated  until  the  objective  function  is  minimized  to  a 
de  s i red  to 1 e  ranee . 


2.  GENERAL  DESCRIPTION  OF  INPUT 

The  following  sections  present  a  general  description  of  the  different 
types  of  data  required  by  the  analysis  and  optimization  programs,  the 
function  of  specific  variables,  and  the  format  of  different  types  of  data. 

In  general,  the  data  has  been  broken  into  separate  blocks  according  to 
function.  These  blocks  are  then  combined  in  the  sequence  shown  in  Figure  2 
to  form  a  complete  input  deck.  Each  of  these  blocks  will  be  described 
separately.  In  some  instances  these  blocks  are  broken  Into  smaller  sub¬ 
blocks  which  are  discussed  separately.  A  detailed  definition  of  each 
variable  is  given  in  the  next  chapter.  Sample  data  sets  are  given  in 
Appendix  A. 
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a.  Executive  Control  Flags 

Two  flags,  FESO  and  FFPE,  are  used  to  select  which  analysis  code  will 
be  used.  The  flag,  OPTM,  is  used  to  select  the  optimization  code.  F2DBLS 
is  set  when  data  for  use  in  the  TWODBL  boundary  layer  program  is  to  be  out¬ 
put  . 

b.  Airfoil  Section  Input 

The  section  geometry  is  input  by  the  same  routine  for  both  the  ESD  and 
the  FPE  codes.  The  values  of  X/C  and  Z/C  for  each  spanwise  control  station 
can  be  input  from  cards  or  from  a  mass  storage  file. 

Streamwise  airfoil  sections  may  be  defined  at  up  to  11  arbitrary  span 
stations.  A  linear  variation  of  coordinates  is  assumed  between  the  input 
stations.  The  first  span  station  must  be  the  wing  root  and  the  progression 
must  be  monotonic  outboard.  Input  for  at  least  two  span  stations  is 
required.  The  flag,  AFIN,  controls  whether  the  section  geometry  will  be 
input  from  cards  or  mass  storage. 

Up  to  90  ordinates  each  may  be  used  to  define  the  upper  and  lower  sur¬ 
face  of  an  airfoil.  Input  coordinates  must  be  normalized  by  the  local  chord 
(X/C, Z/C).  In  the  ESD  code,  computed  airfoil  surface  slopes  may  be  smoothed 
KSMTHS  times  if  required. 

When  inputting  the  airfoil  ordinates  at  each  span  station  from  cards, 
if  the  airfoil  section  is  the  same  as  that  at  the  previous  span  station,  a 
flag,  I  SAME,  is  set  equal  to  T  and  the  program  automatically  uses  the 
previously  defined  airfoil.  If  the  section  at  the  current  span  station 
is  different  from  that  at  the  previous  station  ISAME  is  set  to  F  and  new 
ordinates  are  input.  This  procedure  is  followed  at  each  subsequent  span 
station.  If  the  airfoil  sections  are  input  from  a  storage  file,  the 
ordinates  at  every  span  station  must  be  defined  explicitly. 
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c.  ESO  Input 

The  punched  card  input  data  for  the  ESD  code  are  described  in  this 
section.  Definitions  of  the  input  variables  and  the  required  input  format 
are  presented  in  Subsection  3- 

Program  input  is  divided  into  five  primary  blocks:  program  control, 
wing  planform  geometry,  mesh  skewing  control  parameters,  skewed  mesh 
generation  parameters,  and  Cartesian  mesh  parameters.  These  blocks  are 
delineated  in  Subsection  3  and  briefly  described  below. 

Program  Control 

The  control  section  primarily  governs  the  selection  of  various  program 
options.  Initial  conditions,  required  by  the  solution  technique,  are  defined 
by  either  zeroing  the  potential  function  and  circulation  (IDISK  =  F),  or  by 
reading  these  data  from  a  previous  solution  stored  on  unit  10  (IDISK  *  T) . 

The  latter  option  can  greatly  reduce  computation  time.  If  initial  conditions 
are  to  be  obtained  by  interpolating  from  a  coarser  mesh  solution  stored  on 
unit  10,  the  user  sets  MSHINT  equal  to  T.  Solution  data  for  future  restart 
can  be  saved  on  unit  11  by  setting  ISAVE  equal  to  T. 

The  program  allows  the  storage  of  geometric  and  pressure  data  for  machine 
plotting.  However,  machine  plotting  capability  is  not  currently  available 
with  the  present  system. 

Increases  in  computational  efficiency  and  solution  accuracy  can  be  ob¬ 
tained  by  employing  grid  embedding  and  grid  refinement.  The  recommended 
procedure  is  to  compute  a  solution  on  an  initial  coarse,  skewed  (wing 
oriented)  mesh,  interpolate  the  results,  and  resume  computing  on  an  embedded 
grids  system.  The  embedded  grids  system  is  composed  of  a  fine  interior 
skewed  mesh  embedded  within  a  coarse  Cartesian  exterior  mesh. 
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Two  Namelists,  FLAG  and  FLOW,  are  used  to  input  control  parameters 
and  initial  conditions.  Therefore,  only  those  values  that  will  differ  from 
the  specified  default  values  have  to  be  input. 

Wing  Plan form  Geometry 

A  trapezoidal  reference  wing  is  defined  by  parameters  YROOT,  XLER,  XTER, 
YTIP,  XLET,  and  XTET.  The  wing  reference  area,  chord,  and  moment  center  are 
defined  by  SREF,  CREF,  and  XMOM,  respectively.  All  geometric  input  variables 
are  internally  normalized  by  the  root  chord  of  the  trapezoidal  reference  wing 
CROOT  »  XTER  -  XLER. 

The  wing  leading  and  trailing  edges  are  defined  by  (y,x)  pairs  and  end¬ 
point  slopes  (dx/dy).  Up  to  three  curved  or  straight-line  segments  may  be 
used  to  define  the  leading  and  trailing  edges.  The  input  coordinates  of  each 
segment  are  fitted  with  a  cubic  spline  having  imposed  first-derivative 
(dx/dy)  end  conditions.  Only  two  (y,x)  pairs  are  required  to  define  a 
straight-line  segment.  Breaks  in  planform  are  permissible,  but  the  leading 
and  trailing  edges  must  be  single-valued  in  y  and  the  tip  chord  must  be 
finite.  YP  and  THETAP  define  the  spanwise  location  in  fractions  of  the 
span  and  the  twist  angle  of  each  control  section. 

Me3h  Skewing  Control  Parameters 

This  block  of  input  controls  the  spanwise  skewing  of  the  computational 
mesh.  Two  curves  are  defined  in  the  wing  reference  plane  about  which  the 
mesh  lines  are  aligned.  The  curves  represent  constant  values  of  the  trans¬ 
formed  coordinate  £  and  are  defined  as  the  XI  ■  0  (£  »  0)  and  XI  ■  1  (£  -  1) 
lines.  Generally,  these  curves  should  coincide  with  the  wing  leading  and 
trailing  edges  in  order  to  give  a  desirable  chordwise  distribution  of  grid 
points  on  the  wing. 
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The  XI  *  0  and  XI  -  1  curves  are  described  by  (y,x)  pairs  and  end-point 
slopes  (dx/dy)  in  a  manner  similar  to  the  planform  leading  and  trailing  edge 
description.  The  same  number  of  segments  must  be  used  for  the  XI  ■  0  defini¬ 
tion  as  was  used  for  the  leading  edge  description,  and  similarly  for  XI  *  1 
and  the  trailing  edge  description.  The  XI  lines  must  extend  spanwise  beyond 
the  wing  tip  to  the  edge  of  the  computational  boundary  (ETA^^) ,  and  these 
curves  must  not  intersect. 

Options  are  included  in  the  program  so  that  the  user  can  allow  the  XI  - 
curves  to  coincide  with  all  or  some  segments  of  the  previously  defined 
planform.  Each  previously  defined  segment  of  the  wing  leading  and  trailing 
edges  is  examined  in  turn  to  test  for  acceptability  as  XI  -  lines.  If  the 
leading  or  trailing  edge  segment  is  not  a  satisfactory  definition  for  the  XI  - 
line  segment,  then  (y,x)  coordinates  are  read  in  for  the  segment.  If  the 
outboard  segments  of  the  planform  are  used  for  the  y|  -  line  description,  the 
spline-fit  segments  are  extrapolated  to  the  spanwise  edge  of  the  computational 
space.  If  the  wing  taper  ratio  is  too  small  or  the  tip  too  curved  or  the 
outer  spanwise  boundary  located  too  far  outboard,  the  XI  -  curve  could  Inter¬ 
sect  and  cause  the  program  to  fail. 

Skewed  Mesh  Parameters 

This  portion  of  the  input  applies  to  the  skewed,  wing-oriented  grid 
structure  which  is  either  used  alone  (e.g.,  as  a  single  coarse,  initial  mesh) 
or  in  conjunction  with  a  crude  Cartesian  grid  (embedded  grid  structure).  As 
noted  in  subsection  3  if  a  solution  is  to  be  computedon  an  initial  coarse, 

skewed  mesh  and  the  results  interpolated  for  solution  on  the  embedded  grids 
system  (REMESH  *  T),  then  parameters  For  the  initial  coarse  single  mesh  are 
input  first  followed  by  parameters  for  the  fine  mesh  of  the  embedded  grid 
system.  The  input  parameters  and  input  format  are  essentially  the  same  in 


10 


both  cases.  Namelists  GPARM,  XGRID,  YGRID  and  ZGRIO  are  used  to  input  mesh 
control  parameters. 

The  present  program  has  automatic  grid  generation  capability  in  any  or 
all  of  the  three  coordinate  directions.  If  this  option  is  not  desired  then 
the  appropriate  keys  are  set  to  T  as  indicated  in  Subsection  3*  Up  to  90 
streamwise  mesh  points  (XIN),  30  spanwise  mesh  points  (ETA),  and  20  vertical 
mesh  points  (ZT)  can  be  input.  The  coordinate  system  must  be  consistent 
with  the  planform  coordinate  system  (either  dimensional  or  non-dimensional 
coordinates) . 

If  the  streamwise  grid  is  input,  then  JMAX  values  of  the  mesh  co¬ 
ordinates  along  the  centerline  must  be  specified  by  the  user,  beginning 
upstream  and  proceeding  downstream.  The  mesh  should  be  defined  with  the 
largest  density  of  grid  points  on  the  wing.  The  grid  spacing  should  be 
smallest  near  the  wing  leading  and  trailing  edges  and  should  increase  smoothly 
away  from  the  wing.  The  wing  leading  edge  should  be  situated  midway  between 
grid  points. 

Input  for  automatic  streamwise  grid  generation  is  accomplished  via  name- 
list  XGRID.  All  parameters  in  XGRID  have  preset  (default)  values  which 
have  been  chosen  to  minimize  user  input.  Only  those  parameters  whose  default 
values  are  not  satisfactory  need  be  input.  Two  default  values  for  each  vari¬ 
able  are  shown  in  subsection  3-  The  first  value  applies  to  the  initial  coarse 
mesh  (if  used)  and  the  second  to  the  fine  embedded  mesh.  Parameters  XPLE  and 
XPTE  control  grid  stretching  at  the  leading  and  trailing  edges,  respectively. 
If  a  greater  density  of  grid  points  is  desired  at  the  edges  then  a  smaller 
value  should  be  used.  If  automatic  grid  generation  Is  selected,  JMAX  is  com¬ 
puted  internally  as  NXON  +  NXFWD  +  NXAFT.  NXON,  NXFWD,  and  NXAFT  are  the  spec 
fled  number  of  XI  grid  lines  on,  forward,  and  aft  of  the  wing. 
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If  the  spanwise  grid  is  to  be  read  in,  then  KM AX  values  must  be  speci¬ 
fied  by  the  user  beginning  at  the  wing  root  (ETA  *  0.)  and  proceeding  out¬ 
board.  Again,  the  density  of  grid  points  should  be  greatest  on  the  wing, 
and  the  grid  spacing  smallest  at  the  wing  root  and  tip.  The  tip  should  be 
located  midway  between  spanwise  grid  points.  The  index  of  the  first  grid 
point  outboard  of  the  wing  tip  (KTIP)  must  be  input. 

Some  of  the  above  requirements  are  automatically  satisfied  when  using 
the  spanwise  grid  generation  routine.  Namelist  YGRID  contains  the  user  con¬ 
trolled  variables.  Grid  spacing  at  the  tip  is  controlled  by  YPTIP  and  at 
the  root  by  ALFP,  which  is  the  ratio  of  wing  root-to-tip  stretching  factors. 
For  a  greyer  density  of  grid  points  at  the  root  than  at  the  tip,  input 
ALFP<1.0.  KMAX  is  computed  internally  as  NYON  +  NYOFF.  KTIP  is  computed  as 
NYON  *  1. 

Similarly,  the  vertical  coordinates  (ZT)  at  mesh  points  can  be  input  by 
the  user  or  generated  internally.  When  read  in,  LMAX  values  are  given.  The 
wing  is  positioned  vertically  in  the  computational  space  by  specifying  LWINGU, 
the  index  of  the  first  mesh  point  above  the  wing  reference  plane.  The  grid 
density  should  be  greatest  near  the  wing  reference  plane  and  should  decrease 
with  distance  above  and  beiow  the  wing.  Namelist  ZGRID  contains  the  vertical 
grid  control  parameters.  When  the  vertical  grid  is  generated  internally, 

LMAX  -  ( ( NZ  +  l)/2)  X  2  and  LWINGU  -  LMAX/2  +  1. 

The  parameter  MAXIT  controls  the  total  number  of  relaxation  sweeps  of 
the  grid  system.  When  a  solution  Is  being  sought  for  the  embedded  grids 
system,  MAXIT  refers  to  the  total  number  of  fine  mesh  iterations  for  the 
current  run;  i.e.,  when  the  number  of  iterations  on  the  fine  mesh  equals 
MAXIT,  the  run  terminates  (if  convergence  is  not  obtained  during  the  run). 
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Cartesian  Mesh  Parameters 

Basically,  the  input  procedure  for  the  coarse  Cartesian  mesh  (if  re¬ 
quired  -  EXTMSH  =  T)  is  identical  to  that  for  the  skewed  mesh,  except  that 
now  input  for  controlling  mesh  skewing  is  not  required  (Xi  -  lines).  The 
mesh  coordinates  can  be  input  by  the  user  or  automatic  grid  generation 
routines  can  be  called.  Namelist  names  and  namelist  parameter  names  (and 
their  definitions)  are  identical  to  those  for  the  skewed  grid  except  that 
an  "X"  is  appended  to  the  names.  Namelist  GPARM  is  used  for  the  Cartesian 
mesh  parameters  as  well  as  the  fine  mesh.  However,  only  the  parameters 
appended  with  "X"  must  be  input. 

Parameters  MAXITN  and  MAXITX  control  the  relative  proportion  of  itera¬ 
tions  on  the  fine  and  crude  meshes  of  the  embedded  grids  system.  MAXITN  is 
the  number  of  fine  mesh  iterations  to  be  computed  before  computing  (or  re¬ 
suming  computations)  on  the  coarse  mesh.  Likewise,  MAXITX  is  the  number  of 
coarse  mesh  iterations  to  be  computed  before  resuming  computations  on  the 
fine  mesh.  The  total  number  of  iterations  on  the  fine  mesh  is  given  by 
MAXIT,  which  was  specified  in  the  previous  section, 
d.  FPE  Input 

The  detailed  description  of  the  FPE  input  is  given  in  Subsection  3.  The 
initial  size  of  the  computational  grid  is  specified  by  FNX,  FNY,  and  FNZ.  If 
the  parameter  FHALF  is  not  0,  the  grid  will  be  doubled  after  FIT  iterations. 
Card  5  is  repeated  once  for  each  doubling  of  the  grid.  FCONT  controls 
whether  the  program  computes  an  initial  guess  or  restarts  from  a  stored 
solution.  ZS,  XL,  YL  and  CHD  are  used  to  define  the  shape  of  the  planform. 

XL  and  YL  define  the  coordinates  of  the  leading  edge.  Both  the  leading  and 
trailing  edges  can  have  kinks,  because  the  program  replaces  the  leading  and 
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trailing  edges  as  defined  by  input  with  a  piecewise  straight  line  connecting 
the  nearest  mesh  points  in  the  computational  lattice.  Up  to  1 1  span  stations 
can  be  specified.  Card  13  is  repeated  once  for  each  span  station. 


e.  CONMIN  Input 

The  input  required  when  optimization  is  selected  is  described  in  detail 
in  Subsection  3*  The  various  control  parameters  and  initialization  values 
required  by  CONMIN  are  input  by  the  namelist  OPTIN.  V ( I )  is  the  array  of 
decision  variables  used  to  perturb  the  specified  design  station.  The  IDV 
array  is  used  to  indicate  which  of  the  V  values  are  to  be  used,  i.e.;  for 
IDV  *  1,3,5  -  V ( 1 ) ,  V(3),  and  V{5)  will  be  used.  At  the  present  time,  only 
fifteen  decision  variables  are  used.  CPD  is  the  array  of  desired  pressure 
coefficient  values  that  are  used  to  compute  the  objective  function  and  con¬ 
straint  values.  These  values  should  be  specified  at  x/c  locations  as  close 
as  possible  to  the  x/c  locations  of  the  spanwise  mesh  points  on  the  airfoil 
section  being  optimized. 

Each  of  the  shape  functions  (see  Volume  II,  Part  I,  Page  8)  used  to 
modify  the  base  airfoil  is  multiplied  by  (V(l)  -  1.0)..  Therefore,  initial¬ 
izing  V ( 1 )  ■  1.0  "turns  off"  a  shape  function  until  the  decision  variable, 
V(l),  Is  changed  by  CONMIN. 

In  addition,  various  types  of  constraint  functions  have  been  incor¬ 
porated  into  the  ESO  and  FPE  codes >  Among  these  constraints  are  the  desired 
pressure  coefficients  mentioned  above,  specified  lift  and  moment  coefficients 
and  trailing  edge  thickness.  The  addition  of  other  constraint  functions  will 
require  program  modifications  by  the  user. 


3.  DETAILED  INPUT  VARIABLE  DESCRIPTION 


In  the  following  sections,  the  input  variables  in  each  of  the  five 
blocks  of  data  discussed  in  the  last  chapter  are  defined  on  a  card  by  card 
basis.  Each  block  will  be  treated  separately.  These  blocks  are  separated 
by  title  cards.  In  most  cases,  the  blocks  are  broken  into  sub-blocks  by 
title  cards.  In  the  following  sections,  the  input  cards  are  numbered  in  the 
sequence  they  occur  in  a  block  or  sub-block, 
a.  Executive  Control  Flags 

The  executive  control  flags  are  input  via  the  namelist  CFLGIN.  The 
default  values  are  given  in  parenthesis.  The  control  flag  inputs  are: 


1.  CFLGIN  (NAMELIST  INPUT) 


FESD 

-  T 

selects 

FFPE 

-  T 

selects 

OP  TN 

-  T 

selects 

F2DBLS 

-  T 

signals 

ESD  analysis  code  (F) 

FPE  analysis  code  (F) 

CONMIN  optimization  (F) 

program  to  output  data  to  unit  12  for 


use  in  TWODBL  boundary  layer  program  (F). 


K2DBL 

K0UT8L  (I) 
I  -  1,  20 


■  Number  of  spanwise  stations  at  which  data  for 
TWODBL  will  be  output  (1). 

■  Integer  array  containing  spanwise  indices  of 
stations  at  which  TWODBL  data  is  to  be  written 
(20  *2) 


KOUT  (l)  ■  Integer  array  containing  spanwise  indices  of 

i  *1,  30  stations  at  which  pressure  distribution  data 

will  be  printed/plotted  on  line  printer. 

(30  *  1) 
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b.  Airfoil  Section  Inputs 

Section  ordinates  can  be  input  from  cards  (Unit  5)  or  from  mass  storage 
(Unit  8).  Input  from  mass  storage  assumes  values  of  ordinates  have  been 
written  for  every  spanwise  station.  Mass  storage  input  is  compatible  with 
TV/ODBL  program  output.  The  airfoil  section  inputs  are: 

1.  ATITLE  (1  Card,  8A10) 

80  Character  Airfoil  Section  Description 

2.  AflN  (1  Card,  1 4L5) 

Logical  flag  used  to  indicate  that  section  ordinates  are  to  be 
read  from  cards  (AFIN  *  T)  or  logical  Unit  8  (AFIN  =  F) . 

3.  NPAN ,  I NU , I NL , KSMTHS  (1  Card,  8 1 1 0) 

NPAN  *  No.  of  span  stations  at  which  airfoil  ordinates  are  input 
at  least  two  stations  required  .LE.  11 
INU  ■  No.  of  upper  surface  airfoil  ordinates  .LE.  90 
INL  *  No.  of  lower  surface  airfoil  ordinates  .LE.  90 
KSMTHS  -  No.  of  times  computed  surface  slopes  are  smoothed 

4.  XI NU( I ) , 1  -  I, INU  (INU/8  Cards,  8FI0.0) 

X/C  at  which  airfoil  upper  surface  ordinates  are  input 

5.  XI NL( I ) , I  *  f , I NL( I NL/8  Cards,  8F10.0) 

X/C  at  which  airfoil  lower  surface  ordinates  are  input 

6.  I  SAME  (1  Card,  IH5) 

ISAME  *  F  to  input  airfoil  ordinates 

■  T  use  previously  defined  airfoil 

7.  ZUP(I ,N) , 1-1 , INU  (INU/8  Cards,  8F10.0) 

Upper  surface  airfoil  ordinates/chord  -  required  if  ISAME  *  F 

8.  ZLP( I ,N) , 1*1 , INL  (INL/8  Cards,  8F10.0) 

Lower  surface  airfoil  ordinates/chord  -  required  if  ISAME  *  F 
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Notes : 


A.  I  SAME  must  =  F  for  first  span  station. 

B.  Repeat  cards  6  through  8  NPAN  times  if  AFIN  =  T- 

C.  When  AFIN  -  F,  cards  4,  5,  7  and  8  are  read  from  Unit  8  with 

7  and  8  repeated  NPAN  times.  Card  6  is  omitted. 

c .  ESP  Input 

The  following  input  is  required  when  the  ESP  code  is  selected  (FESD  =  T) . 
The  data  are  broken  into  sub-blocks  in  the  order  that  they  are  input. 

ESP  Control  Flags: 

1.  TITLE  (1  Card.  8A10) 

80  Character  Problem  Title 

2.  FLOW  (Namelist  input,  default  values  are  in  parenthesis) 

MACHNO  *  Free  stream  Mach  number  (.8) 

ALPHAW  *  Wing  reference  plane  angle  of  attack,  degrees  (0.0) 

GAMMA  =»  Specific  heat  ratio  (1.4) 

REACT  =  Regiels  rule  factor  for  modified  slopes  (.01) 

EMEXP(l)  =  Mach  No.  exponent  in  nonlinear  term 
*  2.0  for  Von  Karman-Sprei ter  scaling 
■  1.75  for  modified  Krupp  scaling  (1.75) 

EMEXP(2)  =  Mach  No.  exponent  in  wing  B.  C. 

=  1.0  for  Von  Ka rman-Spre i ter  scaling 
»  -.25  for  modified  Krupp  scaling  (-.25) 

3.  FLAG  (Namelist  input,  default  values  are  in  parenthesis) 

IDISK  =  T  start  from  old  solution  stored  on  Unit  10  (F) 

MSHINT  *  T  initial  conditions  interpolated  from  old  coarser  (F) 
mesh  solution.  IDISK  must  be  true  if  MSHINT  is  true. 
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I  SAVE  *  T  to  save  data  for  future  restart  -  written  on  Unit 

11  (T) 

IPLOT  *  T  to  save  data  for  plotting  -  written  on  Unit  14  (F) 

SOLV  ■  T  for  complete  execution.  Otherwise,  stop  before 

solve  loop  is  entered  (T). 

WBCPRT  -  T  to  output  wing  surface  slopes.  (F) 

FCR  *  T  for  fully  conservative  method  (T). 

ISPAN  **  T  for  inclusion  of  extra  spanwise  terms  in  equation  (T). 

EXTMSH  -  T  to  input  both  fine  interior  mesh  and  coarse  exterior 
mesh  (T) . 

»  F  to  compute  with  fine  interior  mesh  alone 
REMESH  *  T  to  compute  solution  on  initial  mesh,  interpolate,  and 
resume  computing  on  embedded  grids  system 
*  F  to  omit  solution  on  initial  mesh  (T) . 

Planform  Geometry: 

1.  PTITLE  (1  Card,  8A10) 

80  Character  wing  planform  description 

2.  YROOT.XLER.XTER, YTI P.XLET.XTET, SREF,CREF,XM0N  (2  Cards,  6F10.0) 

These  variables  define  trapezoidal  reference  *ing 

YR00T  -  Y  coordinate  of  root 

XIER  ■  X  coordinate  of  L.E.  at  root 

XTER  ■  X  coordinate  of  T.E.  at  root 

YTIP  *  Y  coordinate  of  tip 

XLET  ■  X  coordinate  of  L.E.  at  tip 

XTET  *  X  coordinate  of  T.E.  at  tip 


SREF  ■  Wing  reference  area 


CREF 


Reference  chord 


XMON  *  Moment  reference 

3.  NLES  (1  Card,  8 1 10) 

NLES  -  No.  of  segments  input  to  describe  the  leading  edge 
NLES.LE.3 

4.  NLEI  (1  Card,  8 1 1 0) 

NLEI  *  No.  of  Y,X  pairs  defining  leading  edge  segment  NLEI. LE. 10 

5.  YLE 1(1) ,XLE 1(1), 1*1, NLEI  (NLEI/4  Cards,  8F10.0) 

YLEI.XLEI  ■  Y,X  pairs  defining  the  leading  edge  segment 
At  least  two  pairs  required 

Same  dimensional  system  as  XLER,  etc. 

6.  DXLER.OXLET  (1  Card,  8F10.0) 

DXLER  *  DX/DY  of  L.E.  at  inboard  edge  of  setment 
DXLET  ■  OX/DY  of  L.E.  at  outboard  edge  of  segment 
Note:  Repeat  cards  4  thru  6  NLES  times 

7.  NTES  (1  Card,  8 1 1 0) 

NTES  ■  No.  of  segments  input  to  describe  the  trailing  edge 
NTES.LE.3 

8.  NTEI  (1  Card,  8ll0) 

NTE I  ■  No.  of  Y,X  pairs  defining  trailing  edge  segment 
NTEI. LE. 10 

9.  YTE 1(1), XTE 1(1), 1-1, NTEI  (NTE1/4  Cards,  8F10.0) 

YTEI ,XTEI  ■  Y,X  pairs  defining  the  trailing  edge  segment 
At  least  two  pairs  required 

10.  OXTER, DXTET  (1  Card,  8F10.0) 

OXTER  ■  0X/0Y  of  T.E.  at  inboard  edge  of  segment 
DXTET  ■  0X/0Y  of  T.E.  at  outboard  edge  of  segment 
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NOTE:  Repeat  cards  8  thru  10  NTES  times 
11.  YP(N) ,THETP(N) ,N=1 ,NPAN  (NPAN/4  Cards,  8FI0.0) 

YP  »  Fraction  of  semispan  at  which  airfoils  are  defined 
THETP  *  Twist  Angle,  degrees,  at  YP  -  Positive  is  L.E.  up 
Mesh  Skewing  Parameters: 

1.  STITLE  (1  Card,  8A10) 

80  Character  description  of  skewed  mesh  parameters 

2.  NSO  (i  Card,  8F10.0) 

NSO  *  No.  of  segments  defining  XI  =  0  line 

3.  X0RD{l) ,1=1 ,NS0  (1  Card,  14L5) 

XORD(l)  =  T  Read  in  XI  =0  line  for  Segment  1 

=  F  X I =0  line  for  Segment  1  same  as  wing  L.E.  Seg.  1 
NOTE:  If  XORD(l)  =  F,  skip  cards  4  to  6  for  Segment  1 

4.  NXO  (1  Card,  8 1 1 0) 

NXO  =  No.  of  Y,X  pairs  defining  segment  (NX0.LE.10) 

5.  YXO(I) ,XX0(l) ,  1=1, NXO  (NX0/4  cards,  8F10.0) 

YX0.XX0  =  Y,X  pairs  defining  segment 

6.  DXR0.DXT0  (1  Card,  8F10.0) 

0XR0  =  DX/DY  at  inboard  edge 

DXTO  *  DX/DY  at  outboard  edge 

NOTE:  Repeat  cards  4  thru  6  for  each  X0RD(l)*T 

7.  NTIPLE.NTIPXI  (1  Card,  8ll0) 

Input  NTIPLE.NTIPXI  =  0 

8.  NS  I  ( 1  Card,  8 1 1 0) 

NSI  -  No.  of  segments  defining  XI  =  1  line 
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9.  XIRD(l) ,1=1 , NS  I  (1  Card,  14L5) 


XIRD(l)  =  T  Read  in  Xl  =  l  line  for  Segment  1 

=  F  Xl=l  line  for  Segment  I  same  as  wing  T.E.  Seg.  1 
NOTE:  If  XIRD(l)  =  F,  skip  Cards  10  -  12  for  Segment  1 

10.  NX  I  (1  Card,  8  I  1 0) 

NX  I  =  No.  of  Y,X  pairs  defining  segment  (NX  I . LE .10) 

11.  YXI (I) ,XXI (I) ,1=1 ,NYt  (NX  1/4  Cards,  8F10.0) 

YXI.XXI  =  Y,X  pairs  defining  segment 

12.  OXRI.DXTI  (1  Card,  8F10.0) 

DXRI  =  OX/DY  at  inboard  edge 
DXTI  =  0X/0Y  at  outboard  edge 
NOTE:  Repeat  cards  10-12  for  each  XIRD(l)=T 
Mesh  Parameters  for  Fine  Interior  Mesh  (and  Initial  Coarse  Mesh  if 
REMESH  =  T): 

NOTE:  If  REMESH  =  T,  follow  steps  A  and  B 

If  REMESH  =  F,  do  step  C 

A.  Input  initial  coarse  single  mesh  using  Cards  1-11 
'  B.  Repeat  Cards  1-11  for  fine  mesh 
C.  Input  Cards  I  -  II  for  fine  mesh 

1.  TITLEM  (1  Card,  8A10) 

80  Character  mesh  description 

2.  GPARM  (Namelist  input,  default  values  in  parenthesis) 

MAXIT  ■  No.  of  iterations  for  current  run  (200) 

INCR  *  Iteration  increment  between  write  of  residual  data  (1) 
RSUB  *  Subsonic  relaxation  parameter  (1.8) 

RTEST  ■  Convergence  criterion  of  maximum  potential  error  (.001) 
EPS  =  Coefficient  of  PXT  (I .) 
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KMAX  *  Number  of  spanwise  (ETA)  mesh  points  ( 0) * 

JMAX  ■  Number  of  streamwise  (Xl)  mesh  points  (0)* 

LMAX  *  Number  of  vertical  (ZT)  mesh  points  (0)* 

KTIP  *  First  mesh  point  beyond  wing  tip  (0)* 

LWINGU  »  ZT  mesh  index  of  first  point  above  wing  plane  (0)* 

XMRD  ■  T  read  in  XIN  mesh,  ■  F  compute  XIN  mesh  (F) 

YMRD  *  T  read  in  ETA  MESH,  =  F  compute  ETA  mesh  (F) 

ZMRD  ■  T  read  in  ZT  mesh,  *  F  compute  ZT  mesh  (F) 

NOTE:  If  XMRO  -  T,  skip  Card  3,  If  XMRD  -  F,  skip  Cards  4,5 

3.  XGRiD  (Namelist  input,  default  values  in  (  ),  first  no.  applies  to 
initial  coarse  mesh,  second  no.  to  fine  interior  mesh) 

NXAFT  *  No.  of  streamwise  mesh  points  behind  wing  (5,7) 

NXON  *  No.  of  streamwise  mesh  points  ahead  of  wing  (8/10) 

NXFWD  ■  No.  of  streamwise  mesh  points  ahead  of  wing  (8/1 0) 

ALPXF  *  Oist.  from  L.E.  to  upstream  bdy.  In  root  chords  (3-/. 5) 
ALPXA  *  Dist.  from  T.E.  to  downstream  bdy.  In  root  chords  (4./. 5) 
ALF  *  Location  of  wing  T.E.  between  grid  lines.  If  ALF*0,  TE 
on  grid  line.  ALF*. 5,  midway  between,  etc.  (0./0.) 

XPLE  *  Grid  stretching  factor  at  wing 

XPTE  *  Grid  stretching  factor  at  wing  T.E.  (.2/. 2) 

4.  XIN(j) , J*1 , JMAX  (JMAX/8  Cards,  8F10.0) 

XIN  *  X  mesh  along  centerline 

5.  RX  (1  Card,  8F10.0) 

RX  *  Scale  factor  for  XIN  mesh  (XIN  -  XIN*RX) 

*Values  for  these  parameters  are  computed  internally  when  XMRD,  YMRD,  and 
ZMRD  -  F. 
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NOTE:  If  YMRD  *  T,  skip  Card  6 


If  YMRD  «  F,  skip  Cards  7,8 

6.  YGRID  (Namelist  input,  default  values  in  (  ),  first  no.  applies  to 
initial  coarse  mesh,  second  no.  to  fine  interior  mesh) 

NYON  *  No.  of  spanwise  mesh  points  on  wing  (9/25) 

NYOFF  *  No.  of  spanwise  mesh  points  off  wing  (3/5) 

ALPY  *  Y  dist.  from  wing  tip  to  ETA  bdy.  In  semi-spans  (.5/. 25) 
YPTIP  -  Grid  stretching  factor  at  wing  tip  (.5/. 5) 

ALEP  »■  Ratio  of  wing  root-to-wing  tip  stretch  factor  tl./l.) 

7.  ETA(K) , K=1 , KMAX  (KMAX/8  Cards,  8F10.0) 

ETA  =  Spanwise  coordinate  at  mesh  points 

8.  RY  (1  Card,  8FI0.0) 

RY  *  Scale  factor  for  ETA  mesh  (ETA  *  ETA*RY) 

NOTE:  If  ZMRD  *  T,  skip  Card  9 

If  ZMRD  -  F,  skip  Cards  10,  11. 

9.  ZGRID  (Namelist  input,  default  values  in  (  ),  first  no.  applies 
to  initial  coarse  mesh,  second  no.  to  fine  interior  mesh) 

NZ  ■  No.  of  vertical  mesh  points,  even  (16/20) 

ALPZ  ■  Location  of  outer  ZT  boundary  in  root  chords  (3-/2.) 

ZP  *  Grid  stretching  factor  at  ZT-0.  (.2/. 2) 

10.  ZT(L) ,L-1 ,LMAX  (LMAX/8  Cards,  8F10.0) 

ZT  ■  Vertical  coordinate  at  mesh  points 

11.  RZ  (1  Card,  8F10.0) 

RZ  «  Scale  factor  for  ZT  mesh  (ZT  ■  ZT*RZ) 

Parameters  for  Use  with  Coarse  Exterior  Mesh: 

NOTE:  If  EXTMSH  -  F,  skip  Cards  1  -  11 

(Compute  with  only  one  mesh) 
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,  •  .  li>»  Ununti-i  • 


1.  TITLEM  (1  Card,  8A10) 


80  Character  Coarse  Exterior  Mesh  Description 

2.  GPARM  (Namelist  input,  default  values  in  parenthesis) 

MAXITN  ■  Number  of  iterations  on  fine  mesh  before  computing  on 
outer  mesh  (l) 

MAXITX  *  Number  of  iterations  on  coarse  mesh  before  returning  to 
fine  mesh  (l) 

INCRX  =  Print  increment  for  coarse  mesh  iterations  (1) 

RSUBX  ■  Subsonic  relaxation  parameter  for  coarse  mesh  (1.8) 

EPSEX  =  Coefficient  of  RXT  in  coarse  mesh  calculation  (.5) 

JMAXX  ■  Number  of  streamwise  (XINX)  mesh  points*  (0) 

KMAXX  «  Number  of  spanwise  (ETAX)  mesh  points*  (0) 

LMAXX  »  Number  of  vertical  (ZTX)  mesh  points*  (0) 

KTIPX  =»  First  ETAX  mesh  point  beyond  wing  tip*  (0) 

LWNGUX  ■  ZTX  mesh  index  of  first  point  above  wing  plane  *  (0) 
XMRDX  *  T  read  in  XINX  mesh,  =  F  compute  XINX  mesh  (F) 

YMRDX  ■  T  read  in  ETAX  mesh,  —  F  compute  ETAX  mesh  (F) 

ZMRDX  ■  T  read  in  ZTX  mesh,  =  F  compute  ZTX  mesh  (F) 

NOTE:  If  XMRDX  -  T,  skip  Card  3,  If  XMRDX  *  F,  skip  Cards  4,5 

3.  XGRIDX  (Namelist  input,  default  values  in  parentheses) 


NXONX  - 

(14) 

NXFVDX  » 

See 

(8) 

NXAFTX  - 

Description 

ALPXFX  - 

for 

(5.) 

*  Values  for  these  parameters  are  computed  internally  for  XMRDX,  YMRDX,  and 
ZMRDX  -  F. 
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ALPXAX  = 

Card  3 

(5.) 

ALFX 

Previous 

(0.) 

XPLEX 

Sect i on 

(1.) 

XPTEX 

(.5) 

k.  X I  NX  ( J) , J=1 ,  JMAXX  (JMAXX/8  Cards,  8F10.0) 

XINX  =  Streamwise  coordinate  at  mesh  points 

5.  RX  (1  Card,  8F10.0) 

RX  =  Scale  factor  for  XINX  mesh 
NOTE:  If  YMRDX  =  T,  skip  card  6 

If  YMRDX  =  F,  skip  Cards  7,8 

6.  YGRIDX  (Namelist  input,  default  values  in  parentheses) 


NYONX 

See 

(10) 

NYOFFX  = 

Description 

(5) 

ALPYX 

Card  6 

(1.) 

YPTIPX  = 

Previous 

(.5) 

ALFPX 

Section 

(1.) 

7.  ETAX(K) ,K=1 .KMAXX  (KMAXX/8  Cards,  8F10.0) 

ETAX  =  Spanwise  coordinate  at  mesh  points 

8.  RY  (1  Card,  8F10.0) 

RY  =  Scale  factor  for  ETAX  mesh 
NOTE:  If  ZMRDX  =  T,  skip  Card  9 

If  ZMRDX  =  F,  skip  Cards  10,11 

9.  ZGRIDX  (Namelist  input,  default  values  in  parentheses) 


NZX 

See  Card  9 

(20) 

ALPZX  = 

Previous 

(5.) 

ZPX 

Sect  ion 

(.5) 
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10.  ZTX ( L) ,  L*  1 , LMAXX  (LMAXX/8  Cards,  8F10.0) 

ZTX  -  Vertical  coordinate  at  mesh  points 

11.  RZ  (1  Card,  8F10.0) 

RZ  *  Scale  factor  for  ZTX  mesh, 
d.  FPE  Input 

The  following  data  are  required  by  the  FPE  code.  Separate  cards  or  groups 
of  cards  are  separated  by  a  title  card  which  may  be  blank. 

1.  TITLE  (1  Card,  8A10) 

TEST  CASE  TITLE 

2.  Data  title  card 

3.  FNX,  FNY,  FNZ,  FPLOT,  XSCAL,  PSCAL,  FCONT,  FPRINT  (1  Card,  8F10.0) 

FNX  The  number  of  mesh  cells  in  the  direction  of  the  chord 

used  at  the  start  of  the  calculation. 

FNX  *  0  causes  termination  of  the  program. 

FNY  The  number  of  mesh  cells  in  the  direction  normal  to  the 

chord  and  span. 

FNZ  The  number  of  mesh  cells  in  the  span  direction. 

FPLOT  Controls  generation  of  plots. 

FPL0T»0.  for  a  print  plot  but  no  Calcomp  plot 
at  each  span  station. 

FPLOT*!.  for  both  print  plot  and  a  Calcomp  plot 
at  each  span  station. 

FPL0T-2.  for  a  Calcomp  plot  but  no  print  plot  at 
each  span  station. 

FPL0T-3-  for  a  three  dimensional  Calcomp  plot  only. 

XSCAL,  Control  the  scales  of  the  Calcomp  plots. 

PSCAL 

XSCAL>0.  scales  each  section  plot  to  XSCAL 
XSCAL*0.  scales  each  section  plot  to  5.0 
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FIT 


COVO 


PIO 


P20 


P30 


BETAO 


STRIPO 


FHALF 


The  maximum  number  of  iteration  cycles  which  will  be 
computed . 

The  desired  accuracy.  If  the  maximum  correction  is  less 
than  COVO  the  calculation  terminates  or  proceeds  to  a 
finer  mesh,  otherwise,  the  number  of  cycles  set  by  FIT 
are  completed. 

The  subsonic  relaxation  factor  for  the  velocity  potential. 
It  is  between  1.  and  2.  and  should  be  increased  towards 
2.  as  the  mesh  is  refined. 

The  supersonic  relaxation  factor  for  the  velocity 
potential.  It  is  not  greater  than  1.  and  is  normally 
set  to  1 . 

The  relaxation  factor  for  the  circulation.  It  is 
usually  set  to  I . ,  but  can  be  increased. 

The  damping  parameter  controlling  the  amount  of  added 
4>s  t .  It  is  normally  set  between  0.  and  0.25. 

Determines  the  split  between  horizontal  and  vertical 
line  relaxation  and  is  the  proportion  of  the  total  mesh 
in  which  horizontal  line  relaxation  is  used.  Fastest 
convergence  is  usually  obtained  by  setting  STRIPO  ■  1. 
so  that  horizontal  line  relaxation  is  used  for  the 
entire  mesh.  If  convergence  difficulties  are  encountered 
STRIPO  may  be  reduced  to  some  fraction  between  0.  and  1. 
Determines  whether  the  mesh  will  be  refined.  FHALF*0. : 
the  computation  terminates  after  completing  the  pres¬ 
cribed  number  of  iteration  cycles  or  after  convergence. 
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FHALF^O:  the  mesh  spacing  will  be  halved  after  FIT 
iterations  have  been  run  on  the  crude  mesh  size. 

Card  5  is  repeated  if  FHALF  is  not  0. 

FHAIF<0:  the  interpolated  potential  will  be  smoothed 
(FHALF)  times. 

6.  Data  Title  Card 

7.  FMACH ,  YA.  AL,  CDO,  PWR,  RUN  (1  Card,  8  F10.0) 

FMACH  The  free  stream  Mach  number. 

YA  The  yaw  angle  of  the  wing  in  degrees. 

AL  The  angle  of  attack  in  degrees.  When  the  wing  is  yawed, 

AL  is  measured  in  the  plane  normal  to  the  leading  edge, 
not  in  the  free  stream  direction. 

CDO  The  estimated  parasite  drag  due  to  skin  friction  and 

separation.  It  is  added  to  the  pressure  drag  (sum  of 
vortex  drag  plus  wave  drag)  calculated  by  the  program 
to  give  the  total  drag. 

PWR  Exponent  of  THICK  parameter  used  to  scale  computed  values 

of  XSING  and  YSING. 

RUN  Flag  used  to  control  whether  a  complete  FPE  run  is 

to  be  executed. 

*  0.,  Generate  mesh  only 
■  I . ,  Execute  complete  run 

8.  Data  Title  Card 

9.  ZSYM,  SWEEP1 ,  SWEEP2 ,  SWEEP,  DIHEDI,  DIHED2,  DIHED  (1  card,  8F10.0) 

ZSYM  Determines  whether  to  treat  a  wing  on  a  wall  or  an 
isolated  wing. 

ZSYM«1 . :  the  wing  i s  on  a  wall 
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ZSYM“0.:  the  wing  is  an  isolated  wing  at  a  yaw  angle 
given  by  YAW. 


SWEEP1  Sweep  of  singular  line  at  the  wing  root  if  ZSYM»1 . , 
or  at  the  leading  tip  if  ZSYM»-0. 

SWEEP2  Sweep  of  singular  line  at  the  tip. 

(SWEEP!  and  SWEEP2  are  used  as  end  conditions  for 
spline  fitting  the  x  coordinates  of  the  singular  line.) 

SWEEP  Sweep  of  singular  line  in  the  far  field. 

DIHED1  Oihedral  of  singular  line  at  the  wing  root  if  ZSYM=1 . , 
or  at  the  leading  tip  if  ZSYM=0. 

0IHE02  Dihedral  of  singular  line  at  the  tip.  (DIHED1  and 

DIHED2  are  used  as  end  conditions  for  spline  fitting 
the  y  coordinates  of  the  singular  line.) 

DIHED  Dihedral  of  singular  line  in  the  far  field. 

10.  Data  Title  Card 

11.  SR,  CR,  XR,  BO 2  (1  card,  8F10.0)  •  (Units  must  be  consistent  with  Card  13) 

SR  Reference  wing  area 

CR  Reference  chord  length 

XR  Moment  reference  length 

B02  Semi  span  length 

12.  Data  Title  Card 

13-  ZS,  XL,  YL,  CHD,  THICK,  AL  (NPAN  cards,  8F10.0) 

ZS  Span  location  of  the  section. 

XL,YL  x  and  y  coordinates  of  the  leading  edge. 

CHD  The  local  chord  value  by  which  the  profile  coordinates 

are  scaled. 

THICK  Modifies  the  section  thickness.  The  y  coordinates 
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AL 


are  multiplied  by  THICK. 

The  angle  through  which  the  section  is  rotated  to  intro¬ 


duce  twist.  In  the  case  of  a  yawed  wing,  this  angle  is 
measured  in  the  axis  system  attached  to  the  wing,  not  in 
the  direction  of  the  free  stream. 

Note:  Card  13  is  repeated  once  for  each  span  station, 
e.  CONMIN  Input 

The  following  inputs  are  required  when  optimization  is  selected  (0PTM=T) : 

1 .  Data  Title  Card 

2.  OPTIN  (Namelist  Input,  default  values  are  in  parenthesis) 

NPRINT  (5)  Print  Control.  All  printing  is  done  on  file 
number  6. 

■  0  Print  nothing. 

=  1  Print  initial  and  final  function  information. 

*  2  Print  all  of  above  plus  control  parameters.  Print 
function  value  and  V-vector  at  each  iteration. 

»  3  Print  all  of  above  plus  all  constraint  values, 

numbers  of  active  or  violated  constraints,  direction 
vectors,  move  parameters  and  miscellaneous  inform¬ 
ation.  The  constraint  parameter,  BETA,  printed 
under  this  option  approaches  zero  as  the  optimum 
objective  is  achieved. 

■  M  Print  all  of  the  above  plus  gradients  of  objective 

function,  active  or  violated  constraint  functions 
and  miscellaneous  information. 

■  5  Print  all  of  the  above  plus  each  proposed  design 

during  the  1 -dimens ional  search. 
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NV 


UMAX 


NCON 


NS  IDE 


ICNDIR 


(0)  Number  of  decision  variables,  V(lj,  contained  in 
vector  V. 

(1)  Maximum  number  of  iterations  in  the  minimization 
process.  If  NFDG  =  0,  each  iteration  requires 
one  set  of  gradient  computations  (INFO  =  3  or  b) 
and  approximately  three  function  evaluations  (INFO  = 

1  or  2) .  If  NFDG^O,  each  iteration  requires  approx¬ 
imately  NV  +  3  function  evaluations  (INFO  =  1 
or  2) . 

(0)  Number  of  constraint  functions,  G(J).  NCON  may  be 
zero. 

(1)  Side  constraint  parameter. 

=■  0  The  variables  V(l)  do  not  have  lower  or  upper 
bounds . 

^  0  All  variables  V(l)  have  lower  and  upper  bounds 
defined  by  V L B ( I )  and  VUB(l),  respectively.  If 
any  variables  are  not  bounded  while  others  are,  the 
values  of  the  lower  and  upper  bounds  on  the 
unbounded  variables  must  be  taken  as  very  larqe 
neqative  and  positive  values,  respectively. 

(i.e.  VLB  (I)  =  -1.E+J0,  VUB  (l)  =  l.E+10) 

(NV  +  1)  Conjugate  direction  restart  parameter.  If 
the. function  is  currently  unconstrained,  (all 
G ( J) <CT  or  NCON  and  NSIDE  =  0) ,  Fletcher-Reeves 
conjugate  direction  method  will  be  restarted  with 
a  steepest  descent  direction  every  ICNDIR  iterations. 
If  ICNDIR  =  1 ,  only  steepest  descent  will  be  used. 
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(3)  Scaling  control  parameter.  The  decision  variables 
will  be  scaled  linearly. 

<  0  Scale  variables  V(l)  by  dividing  by  SS  CAL  CO, 
where  vector  SSCAL(l)  is  defined  by  user. 

=  0  Do  not  scale  the  variables. 

>  0  Scale  the  variables  every  NSCAL  iterations. 

Variables  are  normalized  so  that  the  scaled  v(l)  = 

V ( I ) / | V ( 1 )  | .  When  u  sing  this  option  it  is  desirable 
that  NSCAL  =  ICNDIR  if  ICNDIR  is  input  as  nonzero, 
and  NSCAL  =  NV  +  1  if  INCDIR  is  input  as  zero. 

(1)  Gradient  calculation  control  parameter. 

=  0  All  gradient  information  is  provided  by  an  external 
subroutine.  This  information  may  be  calculated 
analytically,  or  by  finite  difference,  at  the  user's 
d i scret ion. 

=  I  All  gradient  information  will  be  calculated  by 

finite  difference  in  CONMIN.  External  subroutine 
provides  only  function  values,  OBJ  and  G(J) ,  J  -  1, 
NCON . 

*  2  Gradient  objective  function  is  provided  by  an  ex¬ 
ternal  subroutine  and  gradients  of  active  and  vio¬ 
lated  constraints  are  calculated  by  finite  dif¬ 
ference  in  CONMIN.  This  option  is  desirable  if  the 
gradient  of  the  objective  function  is  easily  obtained 
in  closed  form,  but  gradients  of  constraint  functions 
G(J),  are  unobtainable.  This  option  may  improve 
efficiency  if  several  variables  are  limited  by  lower 


or  upper  bounds. 


FDCH 


FOCHM 


CT 


CTMIN 


(0.)  Not  used  if  NFDG  ■  0.  Relative  change  in  dec’s  ion 
variable  V ( I )  in  calculating  finite  difference 
gradients.  For  example,  FDCH  ■  .01  corresponds  to 
a  finite  difference  step  of  one  percent  of  the  value 
of  the  decision  variable. 

(0.)  Not  used  if  NFDG  =  0.  Minimum  absolute  step  in 
finite  difference  gradient  calculations.  FDCHM 
applies  to  the  unsealed  variable  values. 

(-.05)  Constraint  thickness  parameter. 

If  CT  <  G(J)  <  | CT [ ,  G(J)  is  defined  as  active. 

If  G ( J )  >  j CT ( ,  G(J)  is  said  to  be  violated.  If 
G(J)>  |CTj  G(J)  is  not  active,  CT  is  sequentially 
reduced  in  magnitude  during  the  optimization  pro¬ 
cess.  If  J  CT |  is  very  small,  one  or  more  con¬ 
straints  may  be  active  on  one  iteration  and 
inactive  on  the  next,  only  to  become  active  again 
on  a  subsequent  iteration.  This  is  referred  to  as 
"zig-zagging"  between  constraints.  A  wide  initial 
value  of  the  constraint  thickness  is  desirable  for 
highly  non-linear  problems  to  reduce  the  zig¬ 
zagging  problem. 

(0.)  Minimum  absolute  value  of  CT  considered  in  the 

optimization  process.  CTMIN  may  be  considered  as 
'numerical  zero',  since  it  may  not  be  meaningful 
to  compare  numbers  smaller  than  CTMIN.  The  value 
of  CTMIN  Is  chosen  to  indicate  that  satisfaction 
of  a  constraint  is  acceptable. 

3* 


CTL 


(-.001)  Constraint  thickness  parameter  for  linear 


CTLMIN 

THETA 


PHI 


■  ■..r.-.aug.irw'u1  <****** 


and  side  constraints.  CTL  is  smaller  than  CT  be¬ 
cause  the  zig-zagging  problem  is  avoided  with 
linear  and  side  constraints. 

(.0005)  Minimum  absolute  value  of  CTL  considered  in 
the  optimization  process. 

(1.)  Mean  value  of  the  push-off  factor  in  the  method  of 
feasible  directions.  THETA  is  called  a  push-off 
factor  because  it  pushes  the  design  away  from  the 
active  constraints  into  th«»  feasible  region.  A 
larger  value  of  THETA  is  desirable  if  the  con¬ 
straints,  G(J),  are  known  to  be  highly  non-linear 
and  a  smaller  value  may  be  used  if  all  the  con¬ 
straints  are  known  to  be  nearly  linear.  The 
actual  value  of  the  push-off  factor  used  in  the 
program  is  a  quadratic  function  of  each  G(j), 
varying  from  0.  for  G(J)  =  CT  to  A.  *  THETA  for 
G(J)  *  | CT | .  A  value  of  THETA  *  0.  is  used  in 
the  program  for  constraints  which  are  identified 
by  the  user  to  be  strictly  linear. 

(0.)  Participation  coefficient,  used  if  a  design  is 
infeasible  (one  of  more  G(J)  >  ) CT | ) .  PHI  is  a 
measure  of  how  hard  the  design  will  be  'pushed* 

1 

towards  a  feasible  region  and  is,  in  effect,  a  j 

i 

penalty  parameter.  If,  in  a  given  problem,  a 
feasible  solution  cannot  be  obtained  with  the 
default  value,  PHI  should  be  increased  and  the 
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DELFUN 


DABFUN 


LINOBJ 


ITRM 


VLB  ( I ) 
1-1,20 


problem  run  again.  If  a  feasible  solution  cannot 
be  obtained  with  PHI  *  100.,  it  is  probable  that  no 
feasible  solution  exists. 

(0.)  Minimum  relative  change  in  the  object  function  to 
indicate  convergence.  If  in  ITRM  consecutive 
iterations,  |l.  -  OBJ(J-l)  /  OBJ ( J) | <  DELFUN  and 
the  current  design  is  feasible  (all  G(J)  <  |CT)  ), 
the  minimization  process  is  terminated.  If  the 
current  design  is  infeasible,  (some  G(j)  >  j CT ) ) , 
five  iterations  are  required  to  terminate  and  this 
situation  indicates  that  a  feasible  design  may  not 
exist. 

(0.)  Same  as  DELFUN  except  comparison  is  on  absolute 

change  in  the  object  function,  | OBJ ( J )  -  OBJ ( J- 1 )  | , 
instead  of  relative  change. 

(0)  Linear  objective  function  identifier. 

=  1  The  objective  (OBJ)  is  specifically  known  to  be  a 

strictly  linear  function  of  the  decision  variables, 
V(l). 

«  0  The  objective  is  a  general  non-linear  function. 

(1)  Number  of  consecutive  iterations  to  indicate  conver 
gence  by  relative  or  absolute  changes  (DELFUN  or 
DABFUN). 

Not  used  if  NSIDE  -  0.  VLB(l)  is  the  lower  allow¬ 
able  value  (lower  bound)  of  variable  V(l).  If  one 
or  more  variables,  V(l),  do  not  have  lower  bounds, 
the  corresponding  VLB(l)  must  be  initialized  to  a 
very  large  negative  number  (i.e.  -1.E+I0). 


VUB(I) 
1  =  1  ,20 


Not  used  if  NSIDE  =  0.  VUB(l)  is  the  maximum 

allowable  value  (upper  bound)  of  variable  V(l).  If 

one  or  more  variables,  V ( 1 ) ,  do  not  have  upper 

bounds,  the  corresponding  VUB(l)  must  be  initialized 

to  a  very  large  positive  number  (i.e.  l.E+10). 

SSCAL(I)  Not  used  if  NSCAL  =  0.  Vector  of  scaling  parameters. 

1=1 ,  20 

The  decision  of  if,  and  how,  the  variables  should 
be  scaled  is  highly  problem  dependent,  and  some 
experimentation  is  desirable  for  any  given  class  of 
problems.  Efficiency  of  the  optimization  process 
can  sometimes  be  improved  if  the  variables  are 
either  normalized  or  scaled  in  such  a  way  that  the 
partial  derivative  of  the  object  function  (OBJ)  with 
respect  to  variable  V(l)  is  of  the  same  order  of 
magnitude  for  all  V(l).  SSCAL(l)  must  be  greater 
than  zero  because  a  negative  value  of  SSCAL(l) 
will  result  in  a  change  of  sign  of  V(l)  and  possibly 
yield  erroneous  optimization  results. 

If  NSCAL>0,  vector  SSCAL  need  not  be  initialized 
since  SSCAL  will  be  defined  in  CONMIN  and  its  asso¬ 
ciated  routines. 

If  NSCAL<0,  vector  SSCAL  is  initialized  in  the 
main  program,  and  the  scaled  variables  V ( I )  * 

V( I )/SSCAL( I ) . 

ISC(I)  Not  used  if  NCON  ■  0.  Linear  constraint  identl- 

1=1,  20 

fication  vector. 
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ALPCON 

BICON 

ALPHAX 

ABOBJ1 

ICONCP 

ICONCL 

ICONCM 

ICNTHK 

CLDSN 

CMDSN 

OBJCON 

3.  Data  Title 

4.  KSD,  KSURF 

KSD 

KSURF 

KCP 


If  constraint  G(J)  is  known  to  be  linear  function 
of  the  decision  variables  V(l),  ISC(j)  should  be 
initialized  to  ISC(J)  *  1.  Identification  of 
linear  constraints  may  improve  efficiency  of  the 
optimization  process  and  is  therefore  desirable, 
but  not  essential. 

If  constraint  G(J)  is  non-linear,  ISC(J)  is  init¬ 
ialized  to  ISC(J)  ■  0.  If  G ( J )  is  not  specifically 
known  to  be  linear,  set  ISC(J)  ■  0. 

(.02)Controls  initial  step  size  in  1-dimensional  search. 
(4.)  Controls  subsequent  1-dimensional  search  step  size. 
Limit  of  change  in  decision  variable. 

Limit  of  change  in  OBJ  function. 

»  1  Selects  Cp  constraints. 

*  1  Selects  wing  Cl  constraint. 

■  1  Selects  wing  CM  constraints. 

■  1  Selects  trailing  edge  thickness  constraint 

Design  C^  for  ICONCL  ■  1. 

Design  C^  for  ICONCM  ■  1. 

(.1)  Scale  factor  for  calculation  of  OBJ  function. 

Card 

KCP,  ILE,  ITE,  NCPl ,  NCP2 ,  NVT  (1  Card,  8110) 

Control  Station  being  designed 

■  1,  upper  surface  is  perturbed 

■  2,  lower  surface  Is  perturbed 

Spanwise  grid  Index  at  which  design  Cp  values  are  com¬ 
pared  with  computed  values 
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I  LE  Initial  chordwise  grid  index  of  range  of  Cp  values 
used  to  compute  ACp  constraints. 

ITE  final  chordwise  grid  index  of  range  of  Cp  values  used 
to  compute  Cp  constraints. 

NCP1  Initial  chordwise  grid  index  of  range  of  Cp  values  used 
to  compute  objective  function. 

NCP2  Final  chordwise  grid  index  of  range  of  Cp  values  used 
to  compute  objective  function. 

NVT  Humber  of  decision  variables  V ( I )  to  be  input. 

Note:  In  the  ESD  code,  the  upper  and  lower  surfaces  have  the 

same  chordwise  indices.  For  the  FPE  code,  the  upper  and 
lower  surfaces  have  different  indices. 

5-  Data  Ti tie  Card 

6.  V ( I )  (NVT/8  cards,  8F10.0) 

Decision  variable  array.  Setting  V ( I )  *  1  will  prevent  a 
decision  variable  from  being  used  until  changed  by  CONMIN 

7 .  Data  Title  Card 

8.  IDV  (NV/8  cards,  8110) 

Decision  variable  index  array.  Use.  to  specify  which  of  the 
V ( I )  values  input  will  be  used  to  turn  on  or  off  the  corres¬ 
ponding  shape  functions. 

9.  Data  Title  Card 

10.  CPD  ( ( ITE- I LE) /8  cards,  8F10.0) 

Array  of  design  Cp  values;  up  to  120  values  can  be  input 

1 1 .  Data  Title  Card 

12.  ZTEMIN,  ZTEMAX  (1  card,  8FI0.0) 

Minimum  and  maximum  trailing  edge  thickness  constraints. 
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4.  DESCRIPTION  OF  OUTPUT 

The  ESD  and  the  FPE  codes  generate  printed  output  and  output  written  to 
mass  storage  or  tape.  The  CALCOMP  plot  routines  from  the  original  FL022 
program  have  been  retained  in  the  FPE  code.  The  ESD  code  writes  pressure 
distribution  data  that  can  be  used  for  post  run  plotting  to  logical  unit. 

Both  the  ESD  and  the  FPE  codes  can  store  data  for  restart  purposes  and  for 
the  TWODBL  boundary  layer  program.  The  following  sections  describe  the  out¬ 
put  that  is  printed  and  stored  by  both  the  ESD  and  the  FPE  codes, 
a .  ESD  Output 

An  abbreviated  example  of  a  typical  ESD  run  with  CONMIN  is  given  in 
Appendix  B.  The  output  from  a  normal  ESD  run  can  be  broken  into  six  major 
sections.  These  sections  are 

Section  I  --  Input  Parameters 

This  section  of  output  is  simply  a  listing  of  all  of  the  input  variables. 

Section  II  --  Computational  Grid 

Information  concerning  the  computational  grid  structure  (skewed  or 
Cartesian  grid)  are  listed  here.  This  includes  the  distribution  of  grid 
points  in  each  of  the  three  coordinate  directions  and,  for  the  skewed  grid, 
the  upstream  and  downstream  boundary  locations  as  a  function  of  span.  The 
listed  grid  coordinates  are  normalized  by  the  planform  root  chord,  CROOT. 

ETA  and  ZT  are  physical  normalized  coordinates  (y,z)  but  XI  is  the  trans¬ 
formed  computational  coordinate,  E,  along  the  centerline  (ETA**0). 

Section  III  --  Wing  Planform  Data 

This  section  gives  data  describing  the  geometry  of  the  wing  planform  and 
the  location  of  the  wing  within  the  computational  mesh.  Also  presented  are 
the  strearrwise  ordinates  of  the  XI  -  0  and  XI  «  1  skewlines  at  each  ETA  grid 
station. 
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Section  IV  —  Wing  Surface  Geometry 

The  next  section  of  output  gives  the  local  wing  chord  and  twist  angle, 
and  wing  surface  ordinates  and  slopes  at  each  ETA  span  station.  This  section 
of  output  can  be  suppressed  by  setting  WBCPRT  equal  to  F  in  the  input. 

Section  V  -  Relaxation  History 

This  section  presents  a  history  of  the  relaxation  process.  The  magni¬ 
tude  and  location  of  the  largest  change  in  the  potential  function  on  both  the 
fine  and  coarse  grids  are  given.  The  change  in  potential  on  the  downstream 
boundary,  the  lift  due  to  circulation  and  the  number  of  supersonic  points 
are  also  given. 

Section  VI  --  Pressure  and  Force  Results 

Wing  surface  pressures,  the  corresponding  local  Mach  numbers,  and  sonic- 
line  coordinates  are  listed  at  each  span  station.  These  pressure  distributions 
are  integrated  chordwise  using  Simpson's  rule  to  give  section  force  and 
moment  coefficients  and  spanwise  loadings.  Total  wing  force  and  moment  co¬ 
efficients  are  then  computed  by  integrating  the  section  data  spanwise.  These 
are  normalized  by  the  input  reference  wing  area  and  chord.  Also  shown  is 
the  total  lift  coefficient  computed  from  the  wing  circulation  distribution, 
F(y).  The  circulation  lift  is  usually  slightly  larger  than  the  integrated 
lift  because  of  "pressure  leakage"  near  the  wing  leading  edge. 

When  a  solution  is  first  computed  on  an  initial  coarse  mesh  (REMESH*T) , 
the  above  sections  of  output  are  listed  first  for  the  initial  solution,  then 
repeated  for  the  solution  on  the  final  mesh.  Also,  when  the  final  mesh  con¬ 
sists  of  a  fine  skewed  grid  embedded  within  a  crude  Cartesian  outer  grid 
(EXTMSH-T),  output  section  II  and  III  are  presented  for  both  grids  of  the 
embedded  grids  system. 
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When  CONMIN  is  cycled  with  the  ESD  code,  additional  output  is  gener¬ 
ated-  This  output  includes  the  current  values  of  the  objective  function, 
the  design  variable  array,  the  modified  upper  and  lower  surface  ordinates 
of  the  design  station  and  the  twist  angle  at  each  control  station.  The 
force  data  described  in  section  VI  is  printed  after  the  final  optimization 
cycle  is  executed. 

Whenever  IPLOT  -  T,  the  ESD  programs  will  write  data  to  unit  14  for 
use  in  a  user  supplied  plotting  program.  The  following  cards  are  required 
to  read  the  data  from  mass  storage. 

READ  (14)  KTMI ,  (ETA(K) ,  K-I.KTMl).  (JLE(K),  K“1 ,  KTM I ) ,  (JTE(K), 

K-l ,  KTMI) 

D010  K-l,  KTMI 
JL  -  JLE(K) 

JT  -  JTE(K) 

READ  (14)  (CPU(K.J),  J  *  JL,JT) 

READ  (14)  (CPL(K.J),  J  -  JL,JT) 

READ  (14)  (XOC(K.J),  J  -  JL,JT) 

10  CONTINUE 
where 

KTMI  »  No.  of  spanwise  grid  stations  minus  one 
ETA  *  Spanwise  coordinate 
JLE  ■  Chordwise  grid  index  of  leading  edge 

JTE  “  Chordwise  grid  index  of  trailing  edge 

CPU  *  Upper  surface  pressure  coefficient 
CPL  •  Lower  surface  pressure  coefficient 
XOC  ■  X/C  values  at  each  span  station 


b.  FPE  Output 

The  FPE  code  generates  printed  output  that  is  similar  to  the  output 
from  the  ESD  code.  An  example  of  this  output  is  given  in  Appendix  B.  The 
first  block  of  data  that  is  printed  is  the  reference  values  of  the  wing  area, 
semispan,  mean  chord,  sweep  and  dihedral  that  are  input  into  the  program. 

This  block  is  followed  by  the  computed  x  and  y  coordinates  of  the  singular 
line  of  the  square  root  coordinate  transformation  at  each  spanwise  control 
station. 

If  optimization  is  not  selected,  the  next  data  printed  describes  the 
computational  mesh  system  being  used.  The  coordinates  of  the  unfolded 
sections  produced  by  the  square-root  transformation  at  the  root  and  tip  are 
printed  and  plotted.  This  data  is  followed  by  the  normal  and  spanwise  cell 
distributions  in  the  square  root  plane  and  the  coordinates  of  the  singular 
line  at  the  spanwise  grid  stations.  If  mesh  halving  is  being  used,  these 
blocks  of  data  are  printed  for  only  the  final  mesh  system. 

The  next  major  block  of  data  displayed  is  the  iteration  history.  The 
maximum  correction  to  the  velocity  potential,  the  maximum  residual  of  the 
difference  equations,  the  circulation  at  the  wing  centerline  and  the  number 
of  supersonic  points  are  printed  at  the  end  of  each  cycle  for  each  mesh 
system  used  in  the  run. 

When  a  convergence  criterion  has  been  satisfied  or  a  specified  number  of 
iterations  has  been  completed,  the  program  prints  the  section  lift,  drag, 
and  moment  coefficients  at  specified  span  stations.  The  pressure  distri¬ 
bution  is  printed  and  plotted  at  equal  intervals  in  the  mapped  plane.  The 
final  blocks  of  data  present  the  force  and  moment  coefficients  for  the  entire 
wing. 
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The  print  generated  when  optimization  is  selected  is  essentially  the 
same  as  in  the  ESD  code.  The  values  of  the  objective  function,  the  decision 
variables,  and  the  ordinates  of  the  span  station  being  optimized  are  given. 

The  FPE  code  can  also  generate  CALCOMP  piots  of  the  pressure  distri¬ 
bution  at  each  span  station  or  a  plot  of  the  three-diroens ional  pressure  dis¬ 
tribution  over  the  upper  and  lower  surfaces  separately. 
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5.  JOB  CONTROL  REQUIREMENTS 
a.  Job  Control  Cards 

The  programs  described  in  this  user's  guide  have  executed  on  two  CDC 
computers,  the  7600  under  the  SCOPE  operating  system  and  the  CYBER  176  under 
the  NOS  operating  system.  The  job  control  information  required  to  run  the 
programs  on  the  CYBER  176  under  the  NOS  system  is  presented  in  this  section. 

The  FORTRAN  code  that  makes  up  the  ESO,  FPE  and  CONMIN  programs  is 
stored  in  three  separate  permanent  mass  storage  files  in  UDDATE  program 
library  format.  The  user  can  therefore  compile  these  three  separate  blocks 
of  code  into  binary  decks  that  can  be  linked  together  in  the  manner  desired 
by  the  user.  The  following  examples  illustrate  the  various  ways  the  ESD  and 
FPE  codes  can  be  executed  with  or  without  CONMIN. 


Sample  decks 


1.  Compile  And  Execute  ESD  or  FPE  From  Update  File 

JOB  card 
USER  card 
CHARGE  card 
CALL,  FTN176 
ATTACH  (OLDPL-ESDCON) 
or 

ATTACH  (OLDPL-FPECON) 

UPDATE  (F) 

FTN  (l ,LCM-I ,L-0) 

RFL  (EC-440) 

LDSET  (PRESET-ZERO) 

LGO. 

7/8/9  -  END  OF  RECORD  Multi  punch 
7/8/9 

DATA  OECK 

6/7/8/9  -  END  OF  INFORMATION  Multi  punch 


45 


2.  Compile  And  Execute  ESD  or  FPE  With  CONMIN  From  Update  Files 

JOB  card 
USER  card 
CHARGE  card 
CALL,  FTN176 
ATTACH  (OLDPL-ESDCON) 
or 

ATTACH  (OLDPL-FPECON) 

UPDATE  (F) 

FTN  (l ,LCM-I ,L-0,B-LG0) 

RETURN  (OLDPL,  COMPILE) 

ATTACH  (OLDPL-CONMINS) 

UPDATE  (F) 

FTN  (I , LCM-I ,L-0,B-LG02) 

RFL  (EC-440) 

LDSET  (PRESET-ZERO) 

LOAD  (LGO, LG02) 

EXECUTE. 

7/8/9 

7/8/9 

7/8/9 

DATA  DECK 

6/7/8/9 


3.  Execute  ESD  or  FPE  From  Stored  Binary  Decks 

Replace  ATTACH  of  OLDPL,  UPDATE,  and  FTN  in  Ex.  1  with 

ATTACH (LGO-ESDCONB) 
or 

ATTACH (LGO-FPECONB) 


4.  Execute  ESD  or  FPE  With  CONMIN  From  Stored  Binary  Decks 

Replace  ATTACH'S,  UPDATE'S,  and  FTN'S  in  Ex.  2  with 

ATTACH (LGO-ESDCONB) 
or 

ATTACH (LGO-FPECONB) 

ATTACH (LG02-C0NM I NB) 


5.  Restart  ESO  or  FPE  From  Stored  Solution 


Insert 

ATTACH (TAPE  1 0-ESDRST) 
or 

ATTACH (TAPE4-FPERST) 

into  the  runstream  prior  to  execution.  In  this  example,  ESDRST 
and  FPERST  are  data  files  saved  from  previous  runs. 

6.  Store  A  Restart  Solution  On  File 

I nsert 

OEFINE(TAPEll-ESDSAV) 

or 

0EFINE(TAPE7-FPESAV) 

into  runstream  prior  to  execution.  DEFINE  is  equivalent  to  doing 
a  REQUEST  and  CATALOG  under  the  SCOPE  operating  system. 

7.  Store  Data  For  Use  In  TVODBL  Program 

I nsert 

DEFINE (TAPE 12-TWODAT) 
into  the  runstream  prior  to  execution. 

b .  Computer  Resource  Requirements 

The  amount  of  computer  time  required  to  run  either  the  ESD  or  the  FPE 
code  with  or  without  CONMIN  will  depend  on  a  number  of  factors  such  as  the 
number  of  iterations  performed  and  the  size  of  the  computational  grid. 

When  linked  with  CONMIN,  the  ESD  code  requires  about  200000g  words  of  main 
memory  to  compile  and  execute  plus  approximately  4l3000g  words  of  large  core 
memory  (LCM)  for  real  time  data  storage.  The  FPE  code  plus  CONMIN  requires 
about  l63000g  words  for  execution  and  431000g  words  of  LCM  for  data  storage. 


SECTION  I  I  I 


VISCOUS  ANALYSIS  PROCEDURE 
(1)  Uses  Three  Programs  —  TWODBL,  LINKUP,  LINKDN 
1.  GENERAL  DESCRIPTION 

The  viscous  analysis  procedure  involves  using  3  programs:  (TWODBL, 

LINKUP,  LINKDN).  TWODBL  is  a  2-D  integral  strip  boundary  layer  program.6 
TWODBL  is  used  to  generate  boundary  layer  displacement  thickness  (6*) 
information  at  the  wing  control  stations  so  that  the  <5*‘s  can  be  subtracted 
from  a  "fluid"  wing  geometry  to  produce  a  solid  wing  geometry  or  to  add 
6*'s  to  the  solid  wing  geometry  to  produce  a  "fluid"  wing.  Inputs  to 
TWODBL  come  from  two  sources:  card  input,  and  file  input  which  has  been 
previously  cataloged.  The  card  input  provides  parameters,  constants,  and 
flags  which  are  needed  to  run  the  boundary  layer  program.  The  file  input 
contains  the  surface  pressures  which  were  computed  and  saved  by  a  previous 
potential  flow  solution. 

The  LINKUP  program  is  used  to  add  6*'s  to  a  solid  wing  geometry  to 
obtain  a  "fluid"  wing.  The  LINKDN  program  is  used  to  subtract  6*'s  from 
a  "fluid"  wing  to  obtain  the  solid  wing  geometry. 

In  order  to  start  the  design  procedure  using  optimization  with  either 
the  FPE  or  ESO  codes,  the  user  must  first  select  a  starting  hard  wing  geometry. 
This  wing  must  then  be  converted  to  an  equivalent  "fluid"  wing.  To  do  this, 
the  hard  wing  must  first  be  analyzed  using  either  the  ESD  or  FPE  potential  flow 
codes.  This  analysis  run  produces  pressure  distributions  over  the  hard  wing 
surface.  These  pressure  distributions  are  then  used  in  a  TWODBL  run  to 
create  a  displacement  thickness  file.  The  user  then  runs  LINKUP  which  accesses 
files  containing  hard  wing  ordinates  and  6*'s  and  adds  the  5*'s  to  the  hard 
wing  ordinates  to  produce  a  file  containing  the  "fluid"  wing  ordinates.  These 
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"fluid"  ordinates  are  used  to  start  the  optimization  design  procedure. 

The  optimization  design  procedure  using  the  inviscid  potential  flow 
codes,  whether  FPE  or  ESD,  results  in  a  wing  geometry  that  produces  the 
specified  design  pressure  distributions.  This  geometry  is  the  equivalent 
"fluid"  wing  geometry  that  the  potential  flow  sees  (hard  wing  +  viscous 
boundary  layer  displacement  thickness).  In  order  to  obtain  a  hard  wing 
geometry  that  can  now  be  used  to  build  a  model,  or  full  scale  wing,  the 
boundary  layer  displacement  thickness  must  first  be  computed  and  then  sub¬ 
tracted  from  the  "fluid"  wing,  at  the  wing  control  stations.  To  do  this, 
the  user  runs  TWOOBL,  and  creates  a  displacement  thickness  file.  The  user 
then  runs  LINKDN,  which  accesses  files  containing  "fluid"  wing  ordinates  and 
6*'s,  and  subtracts  the  6* ' s  for  the  "fluid"  wing  ordinates  to  produce  a 
file  containing  the  hard  wing  ordinates. 

Because  these  viscous  runs  are  done  only  at  the  beginning  and  at  the  end 
of  the  design  process  and  because  it  is  important  that  displacement  thicknesses 
are  added  or  subtracted  properly,  it  is  suggested  that  the  user  carefully 
check  the  TWODBL  run  output. 

NOTE:  TWOOBL  is  used  to  compute  displacement  thicknesses.  These 

computations  are  done  on  a  flat  plate  strip,  the  length  of  the 
strip  being  the  chord  length  which  must  be  input  to  TWODBL. 

2.  RELATION  OF  BOUNDARY  LAYER  STRIP  LOCATIONS  TO  WING  CONTROL  STATION 

LOCATION 

In  general,  the  potential  flow  analysis  codes  compute  surface  pressures 
at  spanwise  locations  that  are  not  exactly  at  the  same  location  as  the  control 
stations  where  airfoil  ordinates  are  prescribed.  The  pressures  are  computed 
at  stations  determined  by  the  grid  distribution  in  the  potential  flow  codes. 
What  the  user  must  do  is  specify  (as  input  to  the  potential  flow  code)  the 
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index  of  the  computed  pressure  station,  that  is  closest  to  each  of  the  control 
stations.  This  is  done  through  an  integer  array,  KOUTBL.  (See  description 
of  potential  flow  code  input).  If  there  are  5  control  stations,  5  numbers 
should  be  input  into  KOUTBL.  Because  the  spanwise  grid  distributions  are 
fairly  dense,  the  difference  in  locations  of  control  stations  and  pressure 
stations  will  be  small,  and  displacement  thicknesses  (corresponding  to  pressure 
stations)  can  be  applied  directly  to  the  control  stations. 

The  tip  control  station  is,  however,  an  exception.  Since  the  potential  codes 
do  not  in  general  compute  accurate  tip  pressures,  the  following  procedure 
should  be  followed.  The  last  index  input  into  KOUTBL  should  correspond  to  a 
pressure  station  in  the  region  of  85%  to  90%  of  the  span.  TWODBL  will  then 
compute  the  displacement  thickness  at  the  spanwise  location  determined  by  this 
index.  The  linking  programs  LINKDN  and  LINKUP  will  then  use  these  displace¬ 
ment  thicknesses  to  adjust  the  airfoil  ordinates  at  the  tip  control  station. 

3.  TWODBL  OUTPUT  NEEDED  TO  RUN  LINKING  PROGRAMS 

There  is  a  minimum  of  TWODBL  output  information  that  the  user  needs  as 
card  input  to  the  linking  programs.  These  inputs  will  be  described  later, 
in  the  section  on  the  linking  programs.  We  will  concentrate  here  on  inputs 
which  require  the  user  to  make  a  judgment  as  to  the  appropriate  values. 

For  both  the  LINKUP  and  LINKDN  programs,  one  of  the  required  inputs  is 
the  location  of  the  separation  point,  if  any.  On  the  TWODBL  output  page 
whose  heading  is  "PRINCIPAL  B0UN0ARY  LAYER  INFORMATION"  will  be  found  a 
number  of  columns.  The  first  column  is  the  point  index.  The  column  X/C 
is  the  streamwise  point  distribution  at  which  the  airfoil  ordinates  are 
given.  The  column  headed  by  FORMI  gives  the  boundary  layer  incompressible 


form  factor.  The  number  of  points  printed  depends  on  whether  and  where 
separation  has  occurred.  For  each  upper  and  lower  surface  strip,  the  user 
will  have  to  input  the  separation  index.  If  there  is  no  separation,  then 
the  index  corresponding  to  X/C  ■  1  will  be  input.  If  there  is  separation, 
the  last  index  printed  will  be  input  on  the  condition  that  the  value  for 
FORMI  is  less  than  2.8.  if  FORM!  is  greater  than  2.8,  then  the  next 
smaller  index  is  the  one  that  will  be  input  into  the  linking  programs. 

Here's  an  example.  Assume  there  are  33  airfoil  defining  points  with 
X/C  »  1.0  at  I  -  33.  For  each  strip  where  33  points  are  printed,  there  was 
no  separation.  Then  the  index  to  use  is  33-  Assume  that  for  two  of  the 
strips,  separation  occurs  at  the  same  streamwise  location,  and  assume  that 
the  last  index  printed  was  I  -  27.  For  one  of  the  strips  the  value  of 
FORMI  at  I  *  27  was  2.78.  Then  for  this  strip  the  index  is  27.  For  the 
other  strip,  the  last  index  printed  was  also  27,  but  the  value  of  FORMI 
at  I  -  2/  woS  '•.82.  Then  the  index  to  use  for  this  strip  is  the  next  lower 
one,  I  »  26. 

There  is  an  exception  to  this  rule.  The  exception  involves  cases 
where  the  displacement  thickness  is  changing  very  rapidly  at  the  point  of 
separation.  In  this  case  it  is  necessary  to  take  the  index  for  separation 
(indicated  by  FORMI),  one  or  two  numbers  smaller.  This  should  be  done 
because  the  separation  streamlines  in  the  linking  program  are  determined 
by  extrapolation.  Taking  the  separation  index  at  a  point  where  the  dis¬ 
placement  thickness  is  changing  very  rapidly,  would  make  such  a  streamline 
extrapolat ion  unrealistic. 
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4.  DESCRIPTION  OF  CARD  INPUTS  TO  TWODBL 

TWODBL  reads  a  Tape  12  file  which  contains  surface  pressure  and  X  locations 
and  writes  a  Tape  10  file  for  use  by  LINKUP  or  LINKDN.  Each  strip  requires  a 
separate  set  of  inputs.  So  if  a  wing  has  5  control  stations,  there  will  be 
two  strips  at  each  control  station  and  two  sets  of  inputs  for  each  control 
station  corresponding  to  the  upper  and  lower  surface  for  a  total  of  10  input 
sets.  The  inputs  start  at  the  inward  part  of  the  wing  and  work  out  to  the 
tip  in  the  following  sequence:  Inboard  upper  surface,  then  inboard  lower 
surface,  then  the  next  station  upper  surface,  then  the  lower  surface  and  so 
on.  A  detailed  description  of  TWODBL  may  be  found  in  Reference  6. 


For  each  strip,  the  following  must  be  input: 

Card 

Title  Card 


Card  #2  8F10.0 

GAM,  R,  PTZ,  TTZ,  UPMACH 
Card  n  915,  F10.5 

NST,  NVP ,  NTURB ,  KPVM,  KEM,  KSMTH ,  KSPLN ,  KLE,  KATCH ,  CTHET 
Card  #4  8F10.0 

DLAM,  TLAM,  DTURB,  TTURB,  RTRAN 
Card  If 5  1 6 1 5 

KPRE ,  KG RAD,  KSOE,  KLAM,  KMAIN,  KPROF 
Card  #6  215,  F10.0 

I  SAME,  I  ADD,  CHORD 


Options  Added 


KPVM  -  6 


X,  y,  0^,  TwaIj  may  be  input  in  4F10.0  format 

T  i,  need  not  be  input,  it  is  assumed  to  be  T  , 
wall  r  total 
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KPVM 


Read 


This  option  must  be  used  for  operation  with  inputs 

generated  by  ESD  or  FPE  codes. 

I  SAME,  I  ADO,  CHORD 
(215,  F10.0) 


If  (SAME  =  1 ,  X  locations  of  input,  and  X  locations  of  airfoil 
defining  points  will  be  read  from  TAPE  12 
If  I  SAME  *  0,  X  locations  will  not  be  read 

IADD  *  Number  of  X  locations  which  correspond  to  the  airfoil 

defining  X/C  locations 

CHORD  =  Chord  length  at  current  input  station. 

For  the  first  set  of  6  input  cards  (inboard  upper  surface,  set  ISAME  * 
1  to  read  X  locations.  For  subsequent  sets,  set  ISAME  =  0. 


Dictionary  of  TWCHBL  Variables 
GAM  specific-heat  ratio,  y 

R  gas  constant,  R,  (ft) (lbf)/(slug) (°R) ;  J/(kg)(K) 

2 

PTZ  inlet  or  upstream  relative  total  pressure  (station  0) ,  P^,  lbf/ft  ; 

N/m2 

TTZ  inlet  or  upstream  relative  total  temperature  (station  0),  T^,  °R;  K 

UPMACH  inlet  or  upstream  Mach  number  relative  to  surface,  MQ 

NST  integer  number  of  input  stations  ( < 1 00)  along  boundary-layer  surface 

NVP  integer  number  of  points  desired  in  velocity  profile  at  each  station 

NTURB  integer  number  of  station,  if  any,  at  which  user  wishes  turbulent 

boundary  layer  to  begin  (NTURB  is  usually  zero,  allowing  program 
to  calculate  position  of  transition  to  turbulent  boundary  layer. 
NTURB  may  also  be  given  any  value  from  1  to  NST.  If  NTURB  ■  1, 
initial  values  must  be  given  for  DTURB  and  TTURB.  If  NTURB  >  1, 
initial  values  may  or  may  not  be  given.) 
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integer  from  1  to  7  indicating  which  form  of  surface  flow  distribu¬ 
tion  is  given  as  inputs.  Note:  KPVM  must  be  set  to  7  when  TVODBL 
is  used  in  the  viscous  linking  procedure. 


Pressure . 1 

Pree-stream  velocity  .  2 

Free-stream  Mach  number  .  3 

Ratio  of  pressure  to  total  pressure  .  4 


Ratio  of  free-stream  velocity  to  free-stream  critical  velocity  5 
(See  Added  Options,  KPVM  =  6; 7) 

integer  (0  to  l)  indicating  which  of  the  two  allowable  sets  of  units 
are  used  in  input: 

English  (pounds  force,  slugs,  feet,  seconds,  degrees 


Rankine,  and  foot-pounds)  .  0 

Metric  (Newtons,  kilograms,  meters,  seconds,  degrees 

Kelvin,  and  Joules)  .  1 


Integer  (0,  1,2,  .  .  .)  indicating  number  of  times  distri¬ 
bution  of  free-stream  velocity  is  to  be  smoothed  prior  to 
computation  of  surface  gradients 

integer  (0  or  l)  indicating  manner  in  which  surface  gradients 


are  to  be  calculated: 

Weighted-difference  technique . 0 

Spline  curve-fit  technique  .  1 


integer  (0  to  1)  indicating  type  of  initial  condition  existing  at 


station  1: 

Stagnation  point  or  initial  values  given . 0 

Sharp  leading  edge . 1 


RTRAN 


KPRE 


KG  RAD 

KSDE 


integer  (0  or  1)  indicating  whether  laminar-boundary-layer  separation 
(if  encountered)  should  reattach  as  a  turbulent  boundary  layer: 

Separation  and  stop . 0 

Reattach  .  I 

real  variable  used  when  KATCH  =  1,  indicating  ratio  of  momentum 
thickness  after  reattachment  to  momentum  thickness  at  laminar 
separation 

initial  displacement  thickness,  if  any,  of  laminar  boundary  layer 
at  station  1,  ft;  m  (DLAM  may  be  zero  or  have  some  finite  value.) 
initial  momentum  thickness,  if  any,  of  laminar  boundary  layer  at 
station  1,  ft;  m  (TLAM  may  be  zero  or  have  some  finite  value.) 
initial  displacement  thickness,  if  any,  of  turbulent  boundary 
layer,  ft;  m  (DTURB  may  be  given  for  station  designated  by  NTURB, 

or  for  station  at  which  transition  is  calculated  by  program.) 
initial  momentum  thickness,  if  any,  of  turbulent  boundary  layer, 
ft;  m  (see  DTURB) 

Momentum  thickness  Reynolds  number  used  to  initiate  computation 
of  the  turbulent  transition  index,  NTURB  (i.e.  NTURB  is  computed 
only  when  Re>RTRAN) . 

integer  (0  or  1)  indicating  whether  printing  of  output  from  PRECAL 
is  desired  (see  OUTPUT): 


Output  suppressed . 0 

Output  printed . 1 


integer  (0  or  1,  see  KPRE)  indicating  whether  printing  of  surface 
velocity  and  Mach  number  is  desired. 

integer  (0  or  1,  see  KPRE)  indicating  whether  printing  of  solutions 
of  laminar  and  turbulent  differential  equations  is  desired. 
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KLAM 


integer  (0  or  1,  see  KPRE)  indicating  whether  printing  of  laminar 
calculations  for  location  of  instability  and  transition  is  desired. 

KMAIN  integer  (0  or  1,  see  KPRE)  indicating  whether  printing  of  principal 
calculated  boundary- layer  parameters  is  desired. 

KPROF  integer  (0  or  1,  see  KPRE)  indicating  whether  printing  of  velocity 
profiles  is  desired. 

X  array  of  X-coord i nates  of  input  stations,  ft;  m 

Y  array  of  Y-coordinates  of  input  stations,  ft;  m 

2  2 

PRES  array  of  static  pressure  P  at  X-Y  input  stations,  Ibf/ft  ;  N/m 

UE  array  of  free-stream  velocities  u  relative  to  surface  at  X-Y 

input  stations,  ft/sec;  m/sec. 

ME  array  of  free-stream  Mach  numbers  M^  relative  to  surface  at  X-Y  input 

stations 

POPTZ  array  of  ratios  of  static  pressure  to  inlet  relative  total  pressure 
P/Pq  at  X-Y  input  stations 

VOVCR  array  of  ratios  of  relative  free-stream  velocities  to  inlet  relative 
critical  velocity  ue/ucr  0  at  X_Y  input  stations  (uc|_  Q  is  the 
speed  of  sound  at  Mach  1,  and  is  only  a  function  of  inlet  relative 
total  temperature.) 


TVAL 


5.  HOW  TO  USE  THE  LINKUP  PROGRAM 

In  order  to  obtain  an  effective  "fluid"  wing  from  a  hard  geometry,  run 
the  LINKUP  program.  The  LINKUP  program  requires  both  card  input,  and  input 
on  a  TAPE7  file.  The  TAPE7  file  is  from  a  previous  TWODBl  run  and  contains 
wing  section  geometries  and  displacement  thickness  nformation.  The  follow¬ 
ing  card  inputs  are  required: 

6.  CARD  INPUT  DESCRIPTION  FOR  LINKUP 
CARD  #1 

NC,  NS  (215) 

NC  =  number  of  X/C  points  at  which  control  station  airfoil  ordinates  were 
specified  (in  the  potential  flow  code). 

NS  =  number  of  surface  strips  which  TWODBL  calculated.  For  example,  for  a 
5  control  station  wing,  with  two  strips  per  control  station,  NS  ■  10. 
CARD  tf2 
NTR  (1615) 

NTR  =  Array  of  indices  of  locations  where  transition  was  specified  in  TWODBL. 
There  will  be  NS  values  input  for  NTR. 

The  sequence  of  indices  is  inboard  control  station  upper 
surface  first,  then  inboard  control  station  lower  surface 
and  so  on,  for  each  control  station  proceeding  out  across 
the  span. 

CARD  #3 
NT  (1615) 

NT  -  Array  of  indices  where  separation  occurs.  See  TWODBL  description  on  how 
to  determine  these  indices-  The  sequence  of  indices  is  the  same  as  for 
NTR. 
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7.  HOW  LINKUP  TREATS  SEPARATION 


A  very  simple  model  is  used  in  LINKUP  to  obtain  "fluid"  ordinates  in 
regions  of  separation.  On  the  upper  surface  the  slope  of  the  fluid  airfoil 
is  taken  to  be  a  constant  in  the  region  aft  of  separation.  On  the  lower 
surface  in  the  cove  region,  the  slope  of  the  fluid  airfoil  is  taken  to  be  a 
constant  aft  of  separation,  until  the  point  at  which  the  difference  between 
the  fluid  and  hard  ordinates  is  equal  to  the  displacement  thickness  at 
separation.  After  this  point,  the  displacement  thickness  is  taken  to  be  a 
constant  and  equal  to  the  displacement  thickness  at  the  separation  point. 

8.  LINKUP  OUTPUT 

LINKUP  produces  no  output  on  unit  6.  All  output  is  written  on  TAPE8. 

This  output  consists  of  control  station  airfoil  geometries  for  the  fluid  wing. 
The  fluid  wing  section  geometries  on  TAPE8  are  in  a  format  suitable  for  input 
into  either  of  the  two  potential  flow  analysis  codes.  TAPE8  should  thus  be 
cataloged  so  that  the  user  doesn't  have  to  punch  the  control  station  fluid 
airfoils.  To  see  LINKUP  output,  copy  TAPE8  to  print  before  cataloging. 

9.  HOW  TO  USE  THE  LINKDN  PROGRAM 

In  order  to  obtain  a  hard  wing  geometry  from  the  fluid  wing  that  comes 
out  of  the  optimization  process  we  run  the  LINKDN  program.  The  LINKDN  program 
requires  both  card  input  and  input  on  a  TAPE7  file.  The  TAPE7  file  is  from  a 
previous  TWODBL  run,  and  contains  wing  section  geometries  and  displacement 
thickness  information.  The  following  card  inputs  are  required  to  run  LINKDN. 

10.  CARD  INPUT  DESCRIPTION  FOR  LINKDN 
CARD  #1 

TOCTE  ( FI  0 . 6) 

TOCTE  is  a  constraint  on  the  trailing  edge  thickness  as  a  fraction  of 

local  chord.  Ouring  the  process  of  obtaining  a  hard  wing  from  the 
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fluid  wing,  the  trailing  edge  thickness  cannot  go  below  TOCTE. 

card  n 

NC,  NS,  LSMTH  (315) 

NC  is  the  number  of  chordwise  X/C  points  in  the  control  station 

airfoil  definition. 

NS  is  the  number  of  surface  strips  which  is  equal  to  twice  the  number 

of  control  stations. 

LSMTH  is  the  number  of  smoothings  to  be  done  on  the  final  hard  wing 
sections . 

CARD  #3 
NT  (1615) 

NT  «  Array  of  indices  of  locations  where  transition  was  specified  in 

TWODBL.  There  will  be  NS  values  input  for  NT.  The  sequence  for 
specifying  NT  is  inboard  control  station  first  -  upper  surface 
index,  then  lower  surface  index;  proceeding  in  this  fashion  inboard 
to  outboard. 

11.  LINKON  OUTPUT 

LINKDN  produces  no  output  on  the  standard  print  file,  UNIT  #6.  All  output 
is  written  on  a  TAPE8  file,  which  should  be  saved  for  subsequent  input  to  the 
potential  flow  codes.  LINKDN  produces  output  identical  in  form  to  that  of 
LINKUP,  which  consists  of  wing  control  station  geometries.  To  see  the  output 
of  LINKDN,  the  user  should  copy  TAPE8  to  a  print  file  before  saving  it. 

12.  JOB  CONTROL  REQUIREMENTS 

TWOOBL  and  the  linking  programs  have  run  on  both  the  CDC  7600  and  the 
CYBER  176  computers.  The  following  sample  decks  illustrate  the  job  control 
cards  required  to  execute  TWODBL  and  the  linking  programs  on  either  computer. 

Sample  Decks: 
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Execute  TW0D8L  from  UPDATE  program  library  and  store  output. 
CYBER  176 


Job  Card 
USER  Card 
CHARGE  Card 
CALL,  FTN 176. 

ATTACH  (TAPE 12  =  TWODINP) 
ATTACH  (OLDPL  =  TWODBL) 
UPDATE,  F. 

FTN,  I,  L=0 ,  LCM= I . 

DEFINE  (TAPE10  =  TWODOUT) 
RFL,  EC= 100. 

LDSET,  PRESET  =  ZERO. 

LGO. 

REWIND,  TAPEIO. 

COPYSBF ,  TAPEIO,  OUTPUT. 
REWIND,  TAPEIO. 

7/8/9  Multipunch 
7/8/9  Mu? t [punch 
INPUT  DECK 
6/7/8/9  Multipunch 


CDC  7600 


Job  Card 

ACCOUNT  Card 

ATTACH  (TAPE12,  TWODINP) 

ATTACH  (OLDPL,  TWODBL) 

REQUEST  (TAPEIO,  *PF) 

UPDATE,  F. 

FTN,  I,  L=0,  LCM-I. 

LDSET,  PRESET  =  ZERO. 

LGO. 

REWIND,  TAPEIO. 

COPYSBF,  TAPEIO,  OUTPUT. 

REWIND,  TAPEIO. 

CATALOG  (TAPEIO,  TWODOUT,  ID=  ) 

7/8/9 

7/8/9 

INPUT  DECK 

6/7/8/9 
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Execute  LINKUP  or  LINKDN  from  UPDATE  Program  Library 
CYBER  176 


Job  Card 
USER  Card 
CHARGE  Card 
CALL,  FTN 176. 

ATTACH  (TAPE7  =■=  TWODOUT) 
DEFINE  (TAPE8  =  LINKOUT) 
ATTACH  (OLDPL  =  LINKS) 
UPDATE,  a. 

FTN,  L=0. 

LDSET,  PRESET  =  ZERO. 
LGO. 

REWIND,  TAPE8. 

COPYSBF ,  TAPE8,  OUTPUT. 
7/8/9  Hultipunch 
*C0MPILE  LINKUP 
or 

*C0MP I LE  LINKDN 
7/8/9 

INPUT  DECK 

6/7/8/9  Multipunch 


CDC  7600 


Job  Card 

ACCOUNT  Card 

ATTACH  (TAPE7,  TWODOUT) 

REQUEST  (TAPE8,  *PF) 

ATTACH  (OLDPL,  LINKS) 

UPDATE,  Q. 

FTN,  L-0. 

LDSET,  PRESET  -  ZERO. 

LGO. 

REWIND,  TAPE8. 

COPYSBF,  TAPE8,  OUTPUT. 

REWIND,  TAPE8. 

CATALOG  (TAPE8,  LINKOUT,  ID  ■*  ) 

7/8/9 

COMPILE  LINKUP 
or 

*C0MPI LE  LINKDN 
7/8/9 

INPUT  DECK 
6/7/8/9 


1 

i 
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SECTION  IV 


WING-PYLON-NACELLE  COMPUTER  CODE 
1.  GENERAL  DESCRIPTION  OF  PROGRAM 

The  transonic  wing-py lon/nacelle  interference  program  (TALA)  consists  of 
three  major  parts.  The  first  part  is  the  ESD  wing  code  which  is  essentially 
the  same  as  the  ESD  wing  code  used  in  the  optimization  design  procedure.  The 
major  differences  in  the  ESD  wing  code  used  in  TALA  are  the  modifications 
made  to  incorporate  the  pylon  and  a  different  input  sequence  for  the  wing  con¬ 
trol  station  ordinates.  The  second  part  of  TALA  are  the  routines  for  solving 
the  nacelle  flow  field.  The  third  part  consists  of  a  driver  program  and 
boundary  condition  transfer  routines,  used  to  cycle  back  and  forth  between 
the  wing-pylon  and  nacelle  solutions.  Fig.  3  shows  the  overall  logic. 

TALA  may  be  operated  in  a  number  of  ways:  I)  A  wing  alone  solution  may 

be  obtained,  2)  A  nacelle  alone  solution  may  be  obtained,  3)  The  nacelle  alone 

solution  may  be  imposed  as  boundary  conditions  on  a  wing  solution.  In  this 
mode,  the  nacelle  boundary  conditions  that  are  imposed  are  not  changed,  and  the 
effect  of  the  wing  on  the  nacelle  is  not  accounted  for.  This  will  be  referred 

to  as  the  non  interacting  mode,  4)  A  mutual  interference  solution  where  the 

nacelle  and  wing  are  updated  with  each  call  to  the  corresponding  program. 

In  this  mode,  both  the  effects  of  the  nacelle  on  the  wing,  and  the  effects  of 
the  wing  on  the  nacelle  are  obtained.  This  will  be  referred  to  as  the  inter¬ 
acting  mode. 

The  nacelle  subroutines  are  set  up  to  solve  a  flow  through  nacelle  geometry 
with  no  center  hubs  or  spinners.  In  the  nacelle  alone  mode,  boundary  conditions 
are  specified  at  upstream  infinity,  downstream  infinity,  and  radial  infinity. 

The  nacelle  solution  is  obtained  on  a  computational  grid  which  maps  the 
infinity  boundary  conditions  to  a  finite  computat ional  domain.  More  detailed 
information  on  grids  and  solution  procedure  can  be  found  in  Ref.  7. 
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2.  GENERAL  DESCRIPTION  OF  CARD  INPUT  SEQUENCE 


Here  we  present  a  short  narrative  description  of  some  of  the  more  important 
input  parameters. 

Namelist  Master  is  Input  First 

This  namelist  controls  the  mode  in  which  the  program  is  run,  and  also 
controls  the  iteration  counts.  If  no  mode  flags  are  specified,  the  default 
mode  will  be  the  interacting  mode. 

NCYCLE  specifies  the  total  number  of  cycles  that  will  be  made  between  the 
nacelle  and  wing-pylon  portions  of  the  program.  From  our  experience  NCYCLE 
should  be  set  to  10  or  larger.  ITERNAC  specifies  the  number  of  nacelle 
solution  iterations  for  each  cycle.  Recommend  ITERNAC  *  30.  ITERATW  specifies 
the  number  of  wing-pylon  solution  iterations  for  each  cycle.  Recommend 
ITERATW  *  30. 

Namelist  Flag  is  Input  Next 

These  inputs  are  used  to  set  up  the  nacelle  computation.  See  detailed 
input  description  for  more  information. 

Namelist  Stretch  is  Input  Next 

These  inputs  are  used  to  set  up  the  nacelle  grid.  See  detailed  input 
description. 

Nacelle  geometry  is  input  next  consisting  of  coordinates  which  either 
have  or  have  not  been  scaled  by  the  nacelle  chord.  There  are  3  sections.  (See 
Fig.  4).  The  coordinates  are  input  8  to  a  card  in  8FI0.0  format.  The  inner 
surface  X  or  X/C  values  are  read  in  first.  Next,  the  inner  surface  section 
ordinates  are  read  in.  The  first  set  of  ordinates  correspond  to  K-l  or  THETA-22. 5°. 
Following  the  first  set  of  ordinates,  the  flag  SAMEIN  is  read.  If  SAMEIN  -  .T. 
then  the  second  set  of  ordinates  is  the  same  as  the  first  set,  and  these 
ordinates  do  not  have  to  be  read.  SAMEIN  is  read  for  each  section  from  K«1  to 
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NTHET.  NTHET  is  set  equal  to  8  inside  the  program.  At  any  K  for  which  the 
section  ordinates  are  not  identical  to  the  preceding  section,  SAHEIN  should 

be  set  to  FALSE  and  the  ordinates  read  in  8F10.0  format. 

After  the  inner  surface  nacelle  geometry  comes  the  input  for  the  outer 
surface  nacelle  geometry  in  a  similar  sequence.  The  outer  surface  can  be 
defined  at  X  or  X/C  values  that  are  different  from  those  used  to  define  the 
inner  surface. 

The  Next  Set  of  Inputs  is  for  the  Wing 

The  wing  inputs  are  essentially  the  same  as  those  described  in  Subsect  ions 
3-b  and  c  of  this  volume  with  the  following  exceptions.  In  the  current  program, 
only  the  fine  grid  system  is  used  and  the  input  card  sequence  for  the  airfoil 
section  ordinates  is  different.  The  airfoil  sections  are  input  after  the 
planform  definition  inputs  and  require  the  following  cards: 

1.  ATITLE  (1  Card,  8A10) 

80  Character  Airfoil  Section  Description 

2.  NPAN,  INU,  INL,  KSMTHS  (1  Card,  8 1 1 0) 

NPAN  -  No.  of  span  stations  at  which  airfoil  ordinates  are  input  - 
at  least  two  stations  required  .LE.  11 

INU  *  No.  of  upper  surface  airfoil  ordinates  .LE.  90 

INL  -  No.  of  lower  surface  airfoil  ordinates. LE.  90 

KSMTHS  =  No.  of  times  computed  surface  slopes  are  smoothed. 

3.  YP(N),  THETP(N),  N-l ,NPAN  (NPAN/4  Cards,  8F10.0) 

YP  *  Fraction  of  semispan  at  which  airfoils  are  defined. 

THETP  ■  Twist  angle,  degrees,  at  YP  -  Positive  IS  L.E.  Up 

4.  XINU(I),  1-1,  INU  (INU/8  Cards,  8F10.0) 

X/C  at  which  airfoil  upper  surface  ordinates  are  input. 
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5.  XINL(I),  1=1,  INL  (INL/8  Cards,  8F10.0) 

X/C  at  which  airfoil  lower  surface  ordinates  are  input. 

6.  ISAME  (1  Card,  1AL5) 

ISAME  =  F  to  input  airfoil  ordinates 

=  T  use  previously  defined  airfoil 

7.  ZUP( I ,N) ,  1=1,  INU  (INU/8  Cards,  8F10.0) 

Upper  surface  airfoil  ordinates/chord  -  Required  if  ISAME  =  F 

8.  ZLP( I ,N) , 1=1 , INL  (INL/8  Cards,  8F1Q.0) 

Lower  surface  airfoil  ordinates/chord  -  Required  if  ISAME  *  F 

NOTES 

A.  ISAME  must  =  F  for  first  YP 

B.  Repeat  Cards  6  through  8  NPAN  times 

The  control  parameter,  AFIN,  is  input  via  namelist  FLAG.  Also,  the  last 
card  in  the  planform  definition  block  has  been  moved  inside  the  airfoil  section 
input  block  shown  above.  After  the  wing  planform  and  section  ordinate  inputs 
is  an  input  flag  PYLN.  If  PYLN  *  T  then  a  pylon  geometry  must  be  input.  If 
PYLN  =  F,  then  no  pylon  computation  is  done,  and  no  pylon  input  geometry  is 
requ i red . 

The  next  card  inputs  ETAP,  and  ALPHP  in  2F10.0  Format.  ETAP  is  the  dimen¬ 
sional  spanwise  distance  at  which  the  pylon  is  located.  ALPHP  is  the  pylon 
toe-in  angle.  ETAP  and  ALPHP  must  always  be  input  even  if  a  wing  alone 
solution  is  being  done  because  ETAP  is  needed  to  set  up  the  grids.  If  NACGEN 
in  namelist  FLAG  was  set  to  true,  it  means  that  a  nacelle-pylon-wing  computational 
grid  will  be  automatically  generated. 

Following  the  cards  containing  the  pylon  geometry  is  a  card  which  inputs 
HNAC  and  XFNAC  in  2F10.0  format.  HNAC  and  XFNAC  give  the  location  of  the  front 
face  of  the  nacelle,  and  must  always  be  input  to  set  up  the  grid.  See  Figure  5. 
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3.  DETAILED  INPUT  DESCRIPTION  FOR  WING-PYLON-NACELLE  CODE 
a .  Control  Flags  and  Initial  Conditions 
NOTE:  Default  values  in  parentheses. 

NAMELIST/MASTER/ 

NCYCLE  ■  #  of  interaction  cycles  between  nacelle  and  wing 
ITERNAC  *  #  of  nacelle  code  internal  solution  iterations  for  each 
cycle  through  nacelle  code 

ITERATW  -  #  of  wing  code  internal  solution  iterations  for  each 


WALONE 


cycle  through  wing  code.  (30) 

Flag  which  determines  whether  a  wing  only  solution  is  wanted(F) 


SAVESOL  *  Flag  which  determines  whether  nacelle  solution  field  will 
be  written  on  tape  file  #9  for  future  restart. 

NALONE  *  Flag  which  determines  whether  a  nacelle  alone  solution 
is  wanted. 

NOINTER  *  Flag  which  determines  if  a  non-interacting  nacelle 

solution  is  wanted  (isolated  nacelle  solution  imposed 
on  wing  as  boundary  conditions  with  no  cycling  between 
nacel 1e  and  wing) . 

NOTE:  If  WALONE  -  .F.,  NALONE  *  .F.,  and  NOINTER  -  .F. ,  then 


the  program  will  run  in  the  interaction  mode. 


NAMELIST/FLAG/ 

ALPHA  *  Nacelle  Angle  of  Attack  (Degrees) 

BETA  *  Nacelle  Yaw  Angle.  Toe-in  is  negative  (Degrees) 

NF  ■  #  of  streamwise  grid  pts.  in  front  of  nacelle 

NA  ■  If  of  streamwise  grid  pts.  on  nacelle 

NB  •  #  of  streamwise  grid  pts.  behind  nacelle 

NE  ■  #  of  radial  grid  pts.  outside  of  nacelle 
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Nl 

ALP 

NALP 

RESTART 

MIN 
WSB 
WSP 
I  TERM 

NESD 

NXBODL 

NXBODU 

CHORD 

NSWCH 

KSMT 

SMBOD 

SMS LOP 

KSMTS 


*  #  of  radial  grid  pts.  inside  nacelle  (14) 

=  Fraction  of  grid  pt.  interval  that  trailing  edge  is  off 

grid  line.  (0.5) 

■  If  NALP  ■  0  then  ALP  is  automatically  set  to  zero 

If  NALP  ■  1  then  ALP  must  be  input  (1) 

■  If  .T.  then  nacelle  solution  is  being  continued  from  a 

previous  solution  (Read  in  on  Tape  File  #8)  (F) 

■  Free  stream  Mach  number  (floating  point  number) 

■  Subsonic  relaxation  factor  (1.6) 

*  Supersonic  relaxation  factor  (0.92) 


*  Nacelle  surface  pressure  information  will  be  printed 
every  I  TERM  iterations. 

*  0  Solve  standard  small  disturbance  equation  (0) 

*  1  Solve  extended  small  disturbance  equation 

*  #  of  nacelle  input  body  pts.  on  inner  surface. 

*  #  of  nacelle  input  body  pts.  on  outer  surface. 

■  Nacelle  chord  length. 

*  F,  program  will  not  switch  to  extended  small  disturbance  equation. 


T,  program  will  switch  to  extended  small  disturbance 
equation  after  ITERESO  iterations.  (F) 

#  of  smoothings  to  be  done  on  nacelle  input  ordinates.  (4) 

.T.  Smooth  nacelle  ordinates  (T) 

.F.  .No  smoothing 

.T.  Smooth  nacelle  surface  slopes  (T) 

.F.  Do  not  smooth 

#  of  smoothings  to  be  done  on  nacelle  surface  slopes.  (4) 
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I  TERES  D  ■  tt  of  iterations  after  which  extended  small  disturbance  terms 

are  turned  on  (if  NSWCH  -  .T.)  (30) 

SCALE  =  .T.  nacelle  geometry  is  input  already 
Scaled  by  the  nacelle  chord. 

*  .F.  input  nacelle  geometry  is  in  physical  units  (con¬ 
sistent  with  units  for  nacelle  chord)  (F) 

NOTE:  Nacelle  geometry  description  must  be  in  units 
consistent  with  the  wing  geometry  description. 


NAMELIST/STRETCH/ 

A1 

*  Axial  grid  stretching  factor 

(0.05) 

A2 

*  Axial  grid  stretching  factor 

(4.0) 

A3 

=*  Axial  grid  stretching  factor 

(0.05) 

A4 

*  Axial  grid  stretching  factor 

(2.5) 

CR 

»  Radial  grid  stretching  factor 

(0.35) 

A3R 

=■>  Radial  grid  stretching  factor 

(0.10) 

A4r 

*  Radial  grid  stretching  factor 

(2.0) 

RS 

*  Nacelle  hi-lite  radius 

(This  is  the  distance  from  the  nace 

1  le  centerl ine  to 

leading  edge  of  the  nacelle.  This 

number  must  be  the 

same  as  the  first  number  input  for 

nacelle  ordinates  at 

X/C  -  0.) 

b. 

Inputs  for  Nacelle  Geometry  Description 

XBOOL  (8F10.0)  Array  of  inner  surface  X  or  X/C  nacelle  defining  points, 
8  to  a  card. 

RBODI  (8F10.0)  Array  of  inner  surface  section  ordinates  at  K*l, 

THETA  ■  22.5°.  These  ordinates  are  defined  with  the 
origin  at  the  nacelle  centerline. 
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The  next  two  cards  control  ordinate  input  for  K=2  to  8  (THETA  from  67.5° 
to  337.5°). 

SAME  IN  (LOGICAL)  =  .T.  The  next  set  of  ordinates  is  the  same  as  the 
previous  set. 

=  .F.  The  next  set  of  ordinates  is  different  from  the 
previous  set,  and  will  be  input. 

RBODI  (8F10.0)  Array  of  inner  surface  section  ordinates  (8  to  a  card) 
at  K>1.  These  must  be  input  for  any  section  where 
SAME  IN  =  .F. 

EXAMPLES:  If  all  eight  nacelle  defining  sections  are  the  same  (nacelle 

is  symmetrical)  then  the  input  will  consist  of  the  initial  section  at  K*l, 
followed  by  seven  cards  inputting  SAMEIN  =  .T.  If  all  eight  nacelle  defining 
sections  are  different,  then  the  input  will  consist  of  the  initial  section  of 
K=1  followed  by  seven  sets  of  SAMEIN  =  .F.,  and  RBODI. 

XBODU  (8F 1 0 . 0 )  Array  of  outer  surface  X  or  X/C  nacelle  defining  points, 

8  to  a  card. 

RBODE  (8F10.0)  Array  of  outer  surface  section  ordinates  at  K«1 , 

THETA  =  22.5°. 

SAMEIN  (LOGICAL)  ■  .T.  The  next  set  of  ordinates  is  the  same  as  the 
previous  set. 

■  .F.  The  next  set  of  ordinates  is  different  from  the 
previous  set,  and  will  be  input. 

RBODE  (8F10.0)  Array  of  outer  surface  ordinates  at  K>1. 

NOTE:  Input  sequence  for  outer  section  ordinates  is  the  same  as  for 
the  inner  ordinates. 
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The  inputs  required  to  compute  the  wing  are  basically  the  same  as  those 
described  previously  for  the  ESD  code.  Only  the  changes  to  the  inputs  to 
accommodate  the  pylon/nacelle  will  be  described  here. 

The  parameter  NACGEN  has  been  added  to  NAMELIST  GPARM. 

NACGEN  »  .T.  The  wing  grid  will  be  automatically  generated  to  accom¬ 
modate  the  pylon/nacelle. 

»  .F.  Standard  ESD  code  grid  procedures  will  be  used. 

NOTE:  It  is  recommended  that  the  code  be  run  with  NACGEN  ■  .T. 

The  following  inputs  have  been  added  immediately  following  the  cards 
containing  the  wing  section  definitions: 

PYLN  (LOGICAL)  *  .T.  Pylon  will  be  computed. 

*  .F.  Pylon  will  not  be  computed,  and  pylon  data  cards 
will  be  omi tted. 

See  Figure  5  for  a  schematic  of  pylon  geometry  setup. 

ETAP,  ALPHP  (2F10.0) 

ETAP  *  Spanwise  pylon  location 

ALPHP  =  Pylon  toe-in  angle  (Degrees;  toe-in  is  positive,  note  dif¬ 
ference  in  sign  from  nacelle  toe-in  angle) 

NSEGL  (MO) 

NSEGL  -  #  of  pylon  leading  edge  segments 

XPL ,  ZPL  (4F10.0) 

XPL  *  X  Coordinate  of  segment  end  point  (Streamwise) 

ZPL  -  Z  Coordinate  of  segment  end  point  (Vertical) 

There  will  be  two  pairs  of  (XPL,  ZPL)  on  each  card  and  one  card  for  each 
segment.  Segment  end  points  will  be  input  starting  at  the  nacelle  and  going  up 
towards  the  wing.  Note:  The  pylon  can  have  at  most  two  leading  edge  segments 
and  two  trail  inn  edge  segments. 
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NSEGT  (110) 


NSEGT  =  #  of  pylon  trailing  edge  segments 
XPT,  ZPT  ( 4 F 1 0 . 0) 

XPT  =  X  Coordinate  of  segment  end  point 
ZPT  =  Z  Coordinate  of  segment  end  point 
There  will  be  two  pairs  of  (XPT, ZPT)  on  each  card,  and  one  card  for  each 
segment.  Segment  end  points  will  be  input  starting  at  the  nacelle  and  going 
up  toward  the  wing. 

Note:  The  zero  reference  point  for  XPL,  XPT,  ZPL,  ZPT  is  the  wing  leading 
edge  at  a  spanwise  position  corresponding  to  the  pylon  location.  Thus,  values 
of  ZPL  and  ZPT  will  be  negative.  XPL,  XPT  will  be  negative  if  they  corres* 
pond  to  a  point  in  front  of  the  leading  edge. 

NOPP  (MO) 

NOPP  *  #  of  points  defining  the  pylon  section. 

XPY  (8F1 0 . 0) 

XPY  =  Array  of  X/C  values  defining  the  pylon  section.  X=0  at  pylon 

leading  edge,  and  C  =  pylon  chord.  There  are  NOPP  numbers,  8  to 
a  card. 

YCOUT  (8F10.0) 

YCOUT  *  Array  of  outboard  pylon  section  ordinates  corresponding  to 
XPY.  There  are  NOPP  numbers,  8  to  a  card. 

YCINB  ( 8 F 1 0 . 0) 

YCINB  =  Array  of  inboard  pylon  section  ordinates  corresponding  to 
XPY.  There  are  NOPP  numbers,  8  to  a  card. 

HNAC,  XFNAC  (2F10.0) 

XFNAC  ■  Distance  from  leading  edge  of  wing  to  face  of  nacelle.  The 

leading  edge  location  is  the  one  corresponding  to  the  spanwise 

71 


, ...  .  tfc*.  >^4. As.  ^  *** »v\.  **  •  * 


location  of  the  pylon.  If  face  of  nacelle  is  in  front  of 
leading  edge,  then  XFNAC  will  be  input  as  a  negative  number. 

HNAC  =  Vertical  distance  from  wing  reference  plane  to  centerline  of 
nacelle.  HNAC  will  be  input  as  a  negative  number. 

NOTE:  Regardless  of  whether  PYLN  is  true  or  false,  the  two  cards  con¬ 
taining  ETAP,  ALPHP  and  HNAC,  XFNAC  must  always  be  input. 

The  remainder  of  inputs  is  the  same  as  previously  described  except  that 
the  exterior  mesh  option  cannot  be  used,  and  if  NACGEN  =  .T.  is  selected,  then 
the  grid  parameters  input  through  NAMELIST  GPARM  will  be  overridden. 
k.  JOB  CONTROL  REQUIREMENTS 
a .  Job  Control  Cards 

The  wing-pylon-nacel le  code  has  been  run  on  two  Control  Data  computers  - 
the  CYBER  176  under  the  NOS  operating  system  and  the  7600  under  the  SCOPE 
operating  system.  Because  of  program  size  restrictions  on  the  7600,  the 
wi ng-py lon/nacel le  code  must  be  segmented  to  run  on  this  computer.  The 
addition  of  segmentation  directives  to  the  7600  runstream  is  the  major  dif¬ 
ference  in  the  job  control  requirements  of  the  two  computers.  The  following 
sample  input  decks  illustrate  the  job  control  cards  required  to  run  the  pro¬ 
gram  on  either  computer. 

Sample  decks: 

1.  Compile  and  execute  from  UPDATE  program  library. 

COC  7600 

Job  Card 

ACCOUNT  Card 

ATTACH  (0LDPL,  TALAF) 

UPDATE,  F. 

FTN,l,L-0,LCM=t,  PL-30000. 

LDSET  (PRESET  -  ZERO) 

SEGL0AD  (B-A) 

LOAD  (LG0) 

N0G0. 
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A. 

EXIT. 

7/8/9  -  Multipunch  End  of  Record 
7/8/9 

SEGMENTATION  DIRECTIVES 

7/8/9 

INPUT  DECK 

6/7/8/9  END  OF  FILE 

CYBER  176 

Job  Card 
USER  Card 
CHARGE  Card 
CALL,  FTN176. 

ATTACH  (OLDPL  =  TALAF) 

UPDATE,  F. 

FTN( I , L=0,  LCM=I,  PL=30000) 

RFL ,  EC=550. 

LDSET  (PRESET  =  ZERO) 

LGO. 

7/8/9 

7/8/9 

INPUT  DECK 

6/7/8/9 


2.  Execute  from  previously  stored  binary  file. 

Replace  ATTACH' s,  UPDATE's  and  FTN's  in  Ex.  1  with 

ATTACH ( LGO, TTB IN)  -  7600 
or 

ATTACH (LGO-TTB IN)  -  CYBER  176 

3-  Restart  from  stored  nacelle  and  wing  solutions  and  store  new  solution 
CDC  7600 

Insert 

REQUEST  (TAPE  1 1 ,  *PF) 

REQUEST  (TAPE9,  *PF) 

ATTACH  (TAPE8,  file,  ID-  );  (Nacelle  Solution) 

ATTACH  (TAPE10,  file,  ID-  );  (Wing  Solution) 

before  execution. 

I nsert 

CATALOG  (TAPE9,  filename,  ID-  );  (Nacelle  Solution) 

CATALOG  (TAPE1 1 ,f i lename,  ID-  );  (Wing  Solution) 

after  execution. 


73 


CYBER  176 


I  nsert 


DEFINE  (TAPE9  =  filename) 
DEFINE  (TAPE 11  =  filename) 
ATTACH  (TAPE8  =  fifename) 
ATTACH  (TAPE1Q  *  filename) 

before  execution. 


Segmentation  Directives 

The  following  segmentation  directives  are  required  to  run  on  the  CDC  7600 


computer. 

TREE  TALA- (NACELLE- (GRID) ,MAI N-( I NPUT) 

GRID  INCLUDE  METRI C .SLOPPY ,SMTH .SPLN2 
NACELLE  INCLUDE  TRIO ,AXI S ,SETI N 
NACELLE  INCLUDE  SHAZ, CUTOUT 
NACELLE  INCLUDE  OUTNAC 

MAIN  I NCLUDE  I NTERP , SOLVE , PYLON , STORE , FARBDY , PF I  NT , WNGBDY . PC  I  NT , 

G  CALC , DBND Y , OUTP , FORCP , FORCE , S I MP , I NPUT 
INPUT  INCLUDE  SETUP.SETUPX.MESH.TCOEF.TCOEFX.MESHIN.WINGCO.PYLCO.SLOPY, 
SPLN1 ,SMTH, IC.GRIDGEN 
NACELLE  GLOBAL  PH  I -SAVE 
TALA  GLOBAL  $Q8. 1 0 . $ , $FCL. C. $ 

GLOBAL  I NTER , I ND I X, METRI CS , CALSUB , SAVEN 
GLOBAL  RPM, RP,RCM,THCM,XPM, TPM, XCM, INDEX 
GLOBAL  FLAGS , RELAXP ,PARM ,W I NG ,GR I DN , LCO 
GLOBAL  LARGN , XYCOE .TEMP , PY I N , ETERMS .JUMP ,W( NGBC 
GLOBAL  $STP. END$ , $  1 0. BUF . $ 

MAIN  GLOBAL  LARGX , EXTER , I NDEXX , LCOX .MESHCO , OLD 
END 


b.  Computer  Resource  Requirements 

The  wing-pylon/nacelle  code  requires  approximately  263000g  words  of  core 
to  load  and  execute  the  program  and  about  505000g  words  of  large  core  memory 
(LCM)  for  run  time  data  storage  on  the  CYBER  176.  Because  of  the  size  of  the 
program,  the  segmentation  described  in  the  preceding  section  must  be  performe 
to  execute  the  program  on  the  CDC  7600  computer.  The  longest  segment  will 
require  about  6 1 OOOg  words  of  small  core  memory  (SCM)  and  about  452000 g  words 
of  LCM. 
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1 


Figure  1  . 


Flow  Chart  of  Design  Procedure 


3.  Wing-Pylon-Nacelle  Solution  Procedure 


WING  LEADING  EDGE 
REFERENCE  POINT 


XFNAC 


APPENDIX  A 


31 


*0?5SS  -25;,1!  *07175  *07465 

‘03636  'till*  *06954  .06529 

.*00232  *  °  *0J30°  *01706 

*0162  _  •003S  "• 0082  .  01872 
-.04630  -.05014  -.05312  -‘o5S73 

-.04744  -.04124  -!o  11% 

-.00226  *C00U  ,0005°  *°6C34 

r 

*0065  .0162  .0242  0319 

*07058  ’°$Ul  *06394  *  . 06590 

*??”!  *??005  *06863  .06625 

1  -03011  *0*353  .01742 


.07465 

.06529 

.01706 

*.01872 


.00034 


.07688 

.06053 

.01201 

-.02530 

-.05625 

-.01883 

-.00005 


.07058 

.03724 

.00178 

.0085 

-.05030 

-.04635 

.00125 

-.00180 


.06625 

.01742 


-.00257 

-.05358 

-.03965 

.00211 


-.01349 

-.05603 

-.03189 

.00198 


-.02385 

-.05763 

-.02362 

.00128 


.01445 

.05966 

.06300 

.01053 


.00000  .00346  .00931  01446 

*06907  ’^68??  ,0557a  *Q59“ 

*03091  #S!5i ?  '°6619  *06300 

:°ooo9  ,0  334  -01659 

.0000  -.0161  -.0239  .  0171 

-*0«n  ?}  *-°*021  ‘  *06049 

*00237  -.02856  -.02095 

-.*00009  ,0°1"'  •°0078  -*00024 

■INC3*  INPUT 

'1ACI!,'0**a**LPH*W.-.66,  « 

n*«T  l,lNG  *3*  PL*NfOPM 

49  838  !*!,,,  S*«12  17.365  1 

49.838  3.0713  5.1316 


0384 

.06758 

.06265 

.01215 

-.03104 

-.05812 

-.01586 

.00024 


.02056 

.06314 

.05855 

.00490 

-.03875 

-.05981 

-.01344 

-.00074 


.07816 

.05520 

.00795 

*.03106 

*.05633 

*.01239 

■,00047 


.04378 

.06878 

.05775 

.00791 

-.03683 

-.05729 

-.00914 

-.00082 


.02708 

.06570 

.05298 

.00152 

•,04472 

*.05788 

*,00733 

’.001S5 


17.365 


,  WBCPRT«. FALSE, 
10.034  11.683 


.07837 

.04937 

.00494 

-.03663 

-.05513 

-.00736 

-.00084 


.04975 

.06967 

.05165 

.00472 

-.04184 

-.05527 

-.00396 

-.00151 


.03389 

.06780 

.04643 

.00032 

-.04947 

*.05459 

-.00262 

■.00117 


.07758 

.04305 

.00300 

-.04169 

-.05223 

-.00371 

-.00098 


•  05450 
.07042 
.04461 
.00255 

*.04633 

•.05172 

>.00061 

•.00178 


.04034 

.06903 

.03889 

,00009 

-.05325 

*.04976 

.00066 

*.00043 


0.0 

.7782 

0.0 

.7782 

4.788 

3.726 

2 

4.7*» 

.5016 

3.726 

.5016 

1 

2 

17.365 

10.034 

0.0 

.3595 

5.441 

.3595 

17.365 

11.683 

3.7642 

-.7235 


3.3166 


2.0964 


INSERT  BLANK  CARD  HERE 


PINE  SKEWED  MESH 
•SPAR*  WAXITM00  I 
•XGRJD  t 
87GRID  | 

(ZGR20  * 
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COARSi  CARTtSIAN  MESH 
IGPARM  S 

AXORIOX  NX0NX«16  ,  MXFWDX*6,  NXAFTXsR,  ALFX«.5,  XPLEXsl.5,  XPTEX«,25 
ALPXAX«2.0  ,  t 
IYGR10X  I 
tZSRIOX  t 
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A. 2  SAMPLE  INPUT  FOR  FPE  +  CONMIN  CODE 


SCFIGIN 

F2DBLS- .T. ,KOUT8L-3,5,9, 15, 16*20,K2DBL  -  4.0PTM-.T., 

FFPE  * . T . .KOUT-3, 5, 9. 15,20,25*20,  S 
CONMIN  NAMELIST 
SOPTIN 

NPRINT-5 ,NV«8, I TMAX-2 , NCON-0 , NS  I OE»0 . IC ND I R-0 , NSCAL «0 , NF DG-0 , 

FOCH- .002 , F  DCHM  = .0002 ,CT*.05,CTL».001 .CTMIN-.01 .CTLMIH* .001 . 

THETA3 1 . , PHI-5. , DEL FUN- .0001 , DAB FUN- 1 . E-06 , L I NOBJ-0 , I TRM-5 . VLB -20* -500. . 
VUB- 20*500. , I  SC -200*0, SSCAL=20*0. , ALPCON- .02 , B 1  CON -4 . 0. CL  DSN- .6. 

ALPHAX- . 1 , 

ABOB01-. 1 , 

CMDSN- - . 1 175.0BJC0N-. 1 . ICONCP-1 , ICONCL-1 . ICONCM-1 , ICNTHK-1 .  I 


KSD 

KSURF 

KCP 

ILE 

ITE 

NCP  1 

NCP2 

NVT 

4 

1 

16 

64 

79 

64 

69 

15 

DESIGN 

VAR IABLES 

1 . 

1  . 

1  . 

I  . 

1  . 

1  . 

1 . 

1  . 

1 . 

1  . 

1  . 

1  . 

1  . 

1  . 

1  . 

DESIGN 

VARIABLES 

USED 

1 

2 

4 

5 

6 

7 

8 

14 

DESIGN 

CP  VALUES 

-.95 

-.95 

-.94 

-.94 

-.94 

-.93 

-.93 

-.92 

-  .92 

-  .91 

-  .90 

-  .89 

-.86 

-  .87 

-.86 

-.85 

2TEMAX 

ZTEMIN 

.014 

.022 

T 

SECTION 

1  ORDINATES 

1 

4 

33 

33 

000000 

.002408 

.009507 

.021530 

.038060 

.059039 

.084265 

.113495 

146447 

. 182803 

.222215 

.264302 

.308658 

.354858 

.402455 

.450991 

. 50000 

.54901 

.59754 

.64514 

.69134 

. 73570 

.77779 

.81720 

.85355 

. 88651 

.91573 

.94096 

.96194 

.  9784’’ 

.99039 

.99759 

1 . 00000 


0. 

.002408 

.009607 

.02153 

.03806 

.059039 

.084265 

.113495 

.146447 

. 182803 

.222215 

.264302 

.308658 

.354858 

. 402455 

.450991 

.5 

.54901 

.59754 

.64514 

.69134 

.7357 

.  7779 

.8172 

.85355 

1  . 

.88651 

.91573 

.94096 

.96194 

. 97847 

.99039 

.99759 

F 

0. 

.01235 

.02066 

.02735 

.03357 

.03978 

.04567 

.05095 

.05547 

.05892 

.06119 

.06232 

.06249 

.06186 

.06061 

.05888 

.05676 

.05428 

.05143 

.04819 

.04454 

.04048 

.0361 

.03158 

.02706 

.0078 

.02273 

.01868 

.01519 

.01256 

.01049 

.009 

.0081 

0. 

-.02261 

-.04028 

-.05571 

-.06835 

-.07846 

- .08606 

-.09176 

-.09577 

-.0982 

-  .0991  7 

-.09894 

-.09694 

- .09324 

- .0883 

-.08192 

- .07473 

- .06688 

-.05873 

-.05032 

-.04225 

-.03537 

-.02913 

-.02353 

-.01875 
-  .00657 

T 

-.01455 

-.01143 

-.009 

-.00758 

-.00696 

-.00655 

- .00653 

F 

0. 

.00818 

.01498 

.02091 

.02649 

.0321 

.03789 

.04363 

.04925 

.05458 

.05947 

.06371 

.06715 

.06965 

.07111 

.07145 

.07071 

.06891 

.0663? 

.06314 

.05956 

.05561 

.05108 

.04583 

.04001 

.01144 

.03388 

.02819 

.02316 

.019 

.01579 

.0135 

.01206 

0. 

- .01324 

-.02459 

-.03379 

-.0408 

-.04653 

-.05162 

-.05622 

-.06028 

-.06358 

-.06601 

-  -  .06744 

- .06777 

- .06697 

- .06489 

-.06135 

-.05607 

-.04913 

- .04106 

-.03266 

-.02473 

-.01787 

-.01252 

- .0089 

-.00678 

-.00546 

C 

-.00538 

-.00407 

-.00303 

-.00271 

-.00321 

- . 004 1 3 

-  .005 

r 

0. 

.00871 

P01587 

.02163 

.02676 

.03153 

.0362 

.04067 

.04498 

.04908 

.05297 

.05651 

.05965 

.06228 

.06421 

.06547 
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.06501 

.06573 

.06466 

.06274 

.05994 

.05623 

.05162 

.04626 

.04045 

.01582 

.03496 

.03005 

.02579 

.02225 

.01949 

.01 752 

.01629 

0. 

-.00812 

-.01506 

-.02068 

-.02546 

-.02948 

- .03314 

-.03646 

-.03944 

-.04208 

-.04421 

-.04572 

- .04646 

-.04603 

- .04444 

-.04146 

- .03721 

-.0316 

-.02484 

-.01814 

- .01 161 

-.00602 

-.00211 

.0006 

.00167 

.00213 

.00203 

.0016 

- . 00002 

-.00214 

- . 00387 

- .00491 

- .00537 


INITIAL 

RUN  IN 

FL022-C0NMIN 

TEST 

F  NX 

FNY 

FNZ 

F  P  LOT 

XSCAL 

PSCAL 

FCONT 

FPRINT 

120. 

20. 

24  . 

0  . 

0. 

0. 

Z. 

0. 

FIT 

covo 

P10 

P20 

P  30 

BETA0 

STR IP0 

F  HAL  F 

100. 

.00004 

1.7 

1 .0 

1  .  4 

0.  1 

1  .0 

0.0 

FVIACH 

YA 

AL 

CDO 

POWER 

RUN 

0.80 

0.0 

0.0 

0.01 

1  .0 

1  .0 

ZSYM 

SWEEP  1 

SWEEP2 

SWEEP 

D1HED1 

DIHED2 

DIHED 

1  .0 

37.93 

26.658 

26.653 

0.0 

0.0 

0.0 

SREF 

CREF 

XREF 

B/2 

1.8971 

.5623 

1.133 

3.374 

ZS/K 

XL 

YL 

CHORD 

THICK 

ALPHA 

0.0 

0.0 

0.0 

1.06123 

1 .0 

4.015 

.3117 

.24288 

0.0 

.  93032 

1  .0 

4.015 

.9302 

.72494 

0.0 

.67043 

1 .0 

3.04874 

3.37392 

1.9518 

0.0 

.3213 

1.0 

1 .02631 

A. 3  SAMPLE  INPUT  FOR  WING-PYLON-NACELLE  CODE 


SMASTER 
SAVESOL-.F. , 

NCYCLE-2, 

ITERNAC- 10. 

ITERATW-t0, 

IPON-.F. . 

SEND 
SF  LAG 

SCALE* . F . , 

MIN-0.80  , 

ITERM* 10. 

US  8* 1 .55. 

WSP-.S2. 

NESD-1  . 

RESTART- .F . , 

ALPHA- -1 . 1 . 

8ETA-0. , 

NF-3  . 

NA-12. 

NB-2  , 

N I  -  7  , 

NE«6. 

NTHET-8, 

ALP-.E  . 

NALP-1  . 

NF-5.NA *24. NB-4.NE-12, 81*14, 
CHORD- 4. $8  , 

NXBODL *17, 

NXBODL - 3 1 , 

NF-7  , 

PI  A  *  2  5  , 

PIB-6, 

USB-1 .5. 

NESD-0 . 

NCFLAT-3 . 

NFLAT-8, 

KSMT-20, 

IFLAT-.F. , 

NSUCH- . F .  , 

KSMTS-4, 

SMSLOP-.T. , 

SEND 

SSTRETCH 

A1-.003, 

A2-2.5  . 

A3- .003 , 

A4-2.E  . 

Al- .05 . 

A2-4  .  , 

A3-.0E  . 

CR-  .  35  . 

A3R-  .  1  . 

A4R-2 .  . 

RS-.7E01  , 

SEND 


0.0 

0.0112 

0.0224 

0.0336 

0.0448 

0.0673 

0.0897 

0.  1345 

0.1793 

0.2242 

0.2446 

0.4484 

1.3450 

2.6900 

4.5800 

0.7801 

0.7472 

0.7336 

0.7235 

0.7154 

0.7026 

0.6928 

0.6792 

0.67H 

0.6673 

0.6669 

0.6669 

0.6669 

0.6669 

0.6669 

T 

T 


0.1121 

4.4840 

0.6852 

0.6669 
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t 

\ 


T 

r 

T 

T 

T 


0.0 

0.0067  0.0211 

0.0379 

0.0674 

0.1010 

0 . 1431 

0.1768 

0.2105 

0.2442  0.2779 

0.3115 

0.3452 

0.3709 

0.4126 

0.4984 

0.5557 

0.6231  0.6904 

0.7578 

0.8252 

0.9064 

1 .2330 

1.6310 

2 .1300 

2.5780  3.0260 

3.3630 

3.8110 

4.2590 

4 . 5800 

0.7801 

0.7954  0.8061 

0.8156 

0.8288 

0.8412 

0.8535 

0.8619 

0.8695 

0 . 876  4  0.8828 

0.8887 

0.8341 

0.8990 

0.9035 

0.9124 

0.9187 

0.9236  0.3272 

0.9294 

0. 9303 

0.9303 

0.9303 

0.9303 

0.9303 

T 

T 

T 

T 

0.9303  0.9303 

0.9144 

0.8469 

0.7467 

0.5869 

1 

r 

T 

T 

ADVANCEO  transonic  technologv 

WING  W3A 

SFLOW 

11ACHNO"  .795,  At.  P  HAW*  - 1  .  1 

,S 

SFIAG 

FCR-.F  .  .EXTMSH-.F.  ,U8CPRT«.F.  .  IDISK-.F.  .AF1N-.T.  ,RF.MF.SH«  .F  .  , 

I  SAVE « 

•  T.  ,S 

WING  PLANFORM 

0.0 

0.0  5.441 

17.366 

10.034 

1 1 .683 

49.838 

3-0713  5.8316 

2 

0.0 

2 

0.0  4.788 

3.726 

.7782 

.  7782 

2 

4.788 

3.726  17.365 

10.034 

.5016 

.5016 

1 

0.0 

2 

~5 . 4  4 1 2  17.365 

11.683 

.3595 

.  3S95 

ATT  FLUID  WING 

5  33  33 

0 

0. 

3.7642  .0924 

3.3166 

.2757 

2.0964 

.6 

.7247 

1  . 

-.7235 

0.00000 

.002410  .00961* 

.02153* 

.038060 

.05904 

.084270 

.113490 

. 146450 

.182800  .22221 

. 26430 

.308660 

. 354860 

.402450 

.45099 

.5000 

.549010  .59755 

.645140 

.69134 

. 7357000 

•>7779 

.81720 

.85355 

.886510  .91573 

.940960 

. 961940 

.978470 

.990390 

.997590 

1  .000 
0.0000 

.002408  .009607 

.021530 

.038060 

.059039 

. 0842G5 

.  113495 

.146447 

.182003  .222215 

.264302 

.308658 

. 354858 

.402455 

.450991 

.  3  000 

.549009  .597545 

.645142 

.691342 

. 735698 

.777785 

.817197 

.353553 

.886505  .915735 

.940961 

.961940 

. 970470 

.990393 

.997592 

l  .000 

F 

.0179 

.0322  .04282 

.05139 

.05839 

.0644. 

.06971 

.07350 

.07702  .07950  .08109 

.05120 

.07996 

.07751 

.07440 

.07066 

.0663*  .06153  .05343 

.05105 

.04544 

.03976 

.03422 

.02891 

.02389  .01913  .01466 

.00066 

.01063 

.00728 

.0045  . 

.00242 

.00  111 

.0179 

-.000)  -.0172 

-.0316 

' .0434 

-.05188 

- .05833 

-.05355 

-.06717  -.06942  -.07011 

- .06953 

- .06743 

-  06403 

- .05948 

-.05380 

-.04739  -.04056  -.03341 

-.02642 

-.02003 

-.01460 

-.01057 

-.00730 

-.00474  -.00296  -.00191 

-.00067 

-.00155 

- .00141 

-.0010.1 

-.00089 

-.0*082 

F 
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AD-A110  232 
UNCLASSIFIED 


LOCKHEED-SEORSIA  CO  marietta  F/G  1/3 

NUMERICAL  AIRCRAFT  DESIGN  USING  3-D  TRANSONIC  ANALYSIS  WITH  OPT— ETC (U) 
AUG  81  R  A  WEED*  A  J  SROKOWSK1  F33615-7B-C-3014 

LG81ER0107-VOL-3-PT-1  AFWAL-TR-B1-3091-VOL-3-PT-  NL 


'I 

1 


.01  3 1  .0306  .0434  .0512 


.07695 

.07937 

.08059 

.08049 

.06625 

.06160 

.05567 

.05122 

.02414 

.01939 

.01455 

.01094 

.00145 

.0131 

.0012 

.0  204 

- 

.03437 

-.06626 

-  .06966 

-.07041 

- 

.06966 

-.04663 

-  .03950 

- .03227 

- 

.02533 

-.00366 

- .00210 

-.00122 

- 

.00087 

-.0*147 

r 

.0162 

.02281 

.02373 

.03649 

.06395 

.06618 

.07175 

.07465 

.07585 

.07314 

.06954 

.06529 

.03633 

.02956 

.02300 

.01706 

.00232 

.0162 

0035 

.0082 

- 

.01872 

-.04630 

- .05014 

-.05312 

- 

.05523 

- .04744 

-.04124 

-.03396 

- 

.02627 

-.00127 

.00011 

.00050 

.00034 

-.00226 

C 

.0085 

0162 

0242 

0319 

.05840 

.06141 

.06384 

.06590 

.07058 

.07005 

.06863 

.06625 

.03724 

.0301 1 

.02353 

.01742 

.00173 

.0085 

-.00257 

-.01349 

- 

.02385 

-.05030 

-.05358 

- .05S03 

- 

.05763 

-.04635 

-  .03965 

- .03189 

- 

.02362 

.00125 

.00211 

.00198 

.00123 

r 

. 00000 

.00346 

.00931 

.01445 

.04620 

.05120 

.05570 

.05966 

.06907 

.06817 

.06619 

.06300 

.03091 

.02334 

.01659 

.01053 

. 00009 

.  0000 

.0161 

.0239 

.0321 

-.05658 

-.05891 

-.06021 

- 

.06049 

-.04351 

-.03626 

-.02356 

- 

.02095 

.00237 

.00188 

.0007  8 

.00024 

- . 00009 
T 


5.358 

0.0 

2 

-3.95 

-.9699 

-2.23 

-  .  4 

-2.23 

-.4 

2 

-.85 

.675 

0.0 

0.6 

2.83 

-.275 

2.83 

-.275 

30 

2.99 

0.0 

0.0 

.0024 

.0095 

.0213 

.  1450 

.1810 

.2200 

.2617 

.  4950 

.5436 

.5916 

.6388 

.9451 

.8777 

.  906  7 

.9524 

0.0 

.006 

.0118 

.0172 

.0373 

.0396 

.0412 

.0423 

.0408 

.0392 

.0372 

.0347 

.0188 

.0154 

.0122 

.0068 

0.0 

-.006 

-.0118 

-.0172 

-.0373 

-.0396 

-.0412 

-.0423 

-.0408 

-.0392 

-.0372 

-.0347 

-.0188 

-  .0154 

- .0122 

-.0068 

0579 

.0637 

.06931 

.07359 

.07918 

. 07698 

.07401 

.07041 

.04566 

.  0400-1 

.02451 

.02919 

.00743 

. 00473 

.06289 

.00181 

- .04464 

-.05253 

- .05901 

- .06316 

- .06749 

-.06395 

-.06913 

- .05332 

- .01905 

-.01362 

- .00932 

-.00603 

- .00091 

-.00109 

-.00128 

-.00142 

.0421 1 

.04799 

.05373 

.05913 

.07688 

.07816 

.07837 

.07758 

.06053 

. 05520 

.04937 

.04305 

.01201 

. 00795 

.00494 

.00300 

-  .02530 

-.03106 

- .03663 

-.04169 

-  .05625 

-.05633 

-.05513 

-  .05223 

- .01883 

- .01239 

-.08736 

-.00371 

-.00005 

- . 0004  7 

-.00084 

-.00098 

0384 

. 04378 

.04975 

.05460 

.06758 

.06878 

.06967 

.07042 

.06265 

.05775 

.05165 

.04461 

.01215 

. 00791 

.00472 

.00255 

- .03104 

-.03683 

-.04184 

- .04633 

- .05B1 2 

-.05729 

- .0557  7 

-.05172 

-.01586 

-.00914 

- .00398 

- .00061 

. 0002  4 

-.00082 

-.00151 

-.00178 

.02056 

.02708 

.03389 

.04034 

.06314 

.06570 

.06780 

.06903 

.05855 

.052911 

.04643 

.03089 

.00490 

.00152 

.00032 

.00009 

- .03875 

-.04472 

-.04947 

-.05325 

- .05981 

- . 05788 

-.05459 

-.04976 

-.01344 

-.00733 

-.00262 

.00066 

- .00074 

-.00155 

-.00117 

-.00043 

0377 

.  3585 

.0933 

.  1124 

3056 

.3513 

.3965 

.  4464 

6845 

.  7284 

.7  700 

.8091 

9877 

1  ,0 

0222 

.0268 

.0309 

.0344 

0428 

.  342*1 

.0427 

.0419 

0320 

.0289 

.0256 

.0223 

0022 

.0006 

.0222 

-.0268 

-.0309 

-.0344 

.0428 

-.0429 

-.0427 

-.0419 

.0320 

-.0289 

-.0256 

-.0223 

.0022 

- .  0006 

88 


1 


r 


-1 .75 


-3.95 

F  F 

a 

i 

F 

FINE  SKEWED 

SGPARM  KMAX-27.NACGEN-.T. , KTIP -23 . YMRD- . F . .ZMRD-. 

MAX  IT* 10  ,0 MAX* 54 , L MAX *20 ,LWINGU*12,S 


.  .XMRD-.F. 


A.  4  SAMPLE  INPUT  FOR  7VODBL  PROGRAM 


ATT  W3 

CL  UPPER 

1 .4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

l 

0  0 

0  1 

0 

i 

33  .4534 

ATT  V3 

CL  LOVER 

1.4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

l 

0  0 

0  1 

0 

0 

33  .4534 

ATT  V3 

ROOT  UPPER 

1  .4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .3837 

ATT  V3 

ROOT  LOWER 

1  .4 

1714.87 

8186. 

498. 

.80 

38 

0  S 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .3837 

ATT  V3 

BREAK  UPPER 

1.4 

1714.87 

8186. 

498  . 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .2847 

ATT  W3 

BREAK  LOWER 

1.4 

1714.87 

3186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .2947 

ATT  W3 

ETA-  .8  UPPER 

1.4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .2235 

ATT  W3 

ETA-  .6  LOWER 

1.4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .2235 

ATT  W3 

ETA- . 85  UPPER 

1  .4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .1666 

ATT  W3 

ETA- . 85  LOWER 

1  .4 

1714.87 

8186. 

498. 

.80 

38 

0  5 

7  0 

2 

0 

0 

0.0 

0.0 

0.0 

0.0 

1 

0  0 

0  1 

0 

0 

33  .1666 

1 . 7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 
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ESD  GRID  AND  GEOMETRY  INFORMATION 


APPENDIX  B 
SAMPLE  OUTPUT 


8.1  SAMPLE  OUTPUT  FOR  ESD  +  CONMIN  CODE 


prno<0»mo«oor«in9>4>OMtr 

■  «<6M*400000000»4I*I«, 
>HOOOOOOOOOOOOOO< 


— ‘OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

♦ooooooooooooooooooooooooooooooooooo 

r^ooooooooooooooooooooooooooooooooooo 

m-OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOQ 


•  N9«4>.4»«f^«««»M««»H««Oir<t«mAO«OOOOOOOOOO0OOOOOOOOOOOOOO 
w,HOo'*>flin«#*n^oo^«^w«»0‘OH«^or*ift«ooooooooooooooooooooooooo 
|.4.<HOOOOOOOOOOOOOOOO«««««4«40«9«OOOOOOOOOOOOOOOOOOOOOOOQO 
100000000000000000000000000*40000000000000000000000000 


oooooooooooooeooooooooooo 


0ifi(MYin«4^M0«0»«ii>M4»»*»N49>MntM0w«inr»4w*ifi®»4Mift*t0p*in000040 

KO»iA<<i(B«*4rA*4mnf(nNNM^<M^«4H»40000ooooooo«««N*4»4<4»<ciMnn(«MMOMi«»r<iioo 

CtOMCIM^^^OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOMMmirtOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


osacM^oo^rKMtMn^oo^Mtaaooooooooooooooooooooooooooooooooooo 

Mon«««0t<fi(Min*<aaMinm«0a^(Mooooooooooooooooooooooooooooooooooo 

Oor«iAna<Ain««Mima«tfta«iMinc(ooooooooooooooooooooooooooooooooooo 

Or*-(M*4000000000000*4(Hr>00000000000000000000000000000000000 

o  ooooooooooooooooooooooooooooooooooo 


or>aoN<|io««»««<n«M«««aa«oAno«r‘«ift« 


00000000000000000000000000 

mooooooooooooooooooooooooo 

sooooooooooooooooooooooooo 

nooooooooooooooooooooooooo 


ooooooooooooooooooooooooo 


IO«0>MO«N400044 

0*400000000000 


00000000000000004 


nW<#OM0«OM»»»»MOa«»4p(KnOOOOOOOOOOOOOOOOOeOOOO 
aopina«o>r^omi»ior>0)««i«)NOAOOOooooooooooooooooooo 
H00^«*4r»0***0t000»*0»f->*4«^00»0000000000000000000000 
•a00««0000044*4l00000«tft0400000000000000000000000 
•  •mMfiet*4«40ooo«*MNfian«oiooooooooeeoooeoeoooooo 


*4000000000000000000000000000004 


04>>4r«a<4raar«a9<M9f«ao«atMHa»>4r>rt«^«aoooooooooooooooooooooooo 

<o«<*iMtciaraa»in«M>nain«f«(<(ninar^a««4««44iMoooooooooooooooooooooooo 

»40*«r*or*ooooor><*«(aoif»ooc>«or*ooooinr-ooMOoooooooooooooooooooooooo 

wo^tfnooo*«oo«n44a*mo4f>«4p-noo4>»v>r*o«>MO«4aiooooooooooooooooooooooooo 

oo*«r«n««raor44Mn«ao<44ia«p>aao«4M(<iaa»oooooooooooooooooooooooo 


>000000000000000000 


h  xoooooom«40oo44C<oooomoo«aoo4<woM004^»*oor*«44rr»o^oo««pia»ooooorav^iftOO 
►  aiMp»^^ooooooooooo-H^^^NMrtMMM»»iaaain«i«a«FkM*Mi4iaaaaa*a»oooo*»pia 


y*3U!{*i**+  ***&&■ 


UPSTREAM  AND  DONNSTREAM  BOUNDARY  LOCATION 


MINS  DESCRIPTION 


T 


477 


9  9  »  0> 


o  o  o  o  o 
o  o  o  o  o 
AnKHAM 


1 


«*t«nw  ooooo 
«  r»  s  o«  ooooo 


Q 

a 


_ lO<  ooooo 

irMomin  «  o  «  «  o 


coco 
m  «a  m 


c  c  c  c  c 
nmnmin 


X 

►  u 

Z  X 


hi  X 
(9  • 


ooooooooooooooooooooooooo 

ooooooooooooooooooooooooo 

ooooooooooooooooooooooooo 

ooooooooooooooooooooooooo 

ooooooooooooooooooooooooo 

ooooooooooooooooooooooooo 


x  x  x  r* 
r-  r*  m  x 
w  «o 


•2  U 
M  % 


hfd 

•2  O  O  < 

OX  0  0  4 

z  •  *  • 

X  X 


OOOOOOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooo 


I-  0«9  •  • 


o 

X  fed 


i  r-r-o  ox  o»  o  ( 


o 

u 

X  *2 
hi  *5 


^•xxr^xr--*©**© 


H 

•a  o  o  oxx  o  9“  X  X  X  m  t 

o©x©c4xor-r-r-«*m©< 


»  r«o«  V( 
I90**C( 


o-^r-xf-xxxxxw^xxxxoxxf^xxoxx 


U  H  •  •  «  I  •  •  •  •  t  •  •  •  . . .  •  •  • 


•2 

X 


o 

K 

X 

X 

X 

0*4 

M 

X 

X 

« 

X 

u 

X 

X 

«« 

X 

X 

o 

o 

M 

X  ♦ 

X 

X 

X  X 

X  u 

X 

u 

X 

u  o 

X 

H 

X 

o 

X 

X 

•4 

X 

• 

u 

w 

X 

X 

M 

o 

X 

H  «4 

«  o 

X 

X 

X 

o  o 

X 

o 

X 

X  o 

X 

u 

X 

x  r» 

94 


2J. 


NXOMX  NXFMDX  NXXFTX  XLFXFX  ALPXXX  XLrX  XPLEX  XPTEX 

1441  5.0000  2.0000  .5000  >.5000  .2500 


COMPUTATIONAL  GRID 


«9 

X 

O 

H 

H 

x 

u 

o 

J 

X 

o 


•  o«A*«o««Aminir»»«o«4m«o«oooooooooo 
^^••^^^^♦♦♦♦•^HnrnsD^ifl^oooooooooo 
xtf)«xox«tfx»or'r*o<*«>a»^oxa>tf)oooooooooo 
t4#*m«oxc««4r»»c«r'ior»«4r'M0o«»«"»fnoooooooooo 
•Jr*iNC*oomr^^«r<r*««o-«r-moor*f«r“Oooooooooo 


f  t 


*4*<nin*rin^coooooooooo 


Otft«rf)rt*«nor«r>««(MCDOOOOOopoooooooo 

K0^»r'^^l/>P«^«9*^P‘^I,*000000000000000 

<orMA«r»a»»r>«Nf»MMooooooooooooooo 

HOOff>^W<m^0i«»0*MO»0O©OOOOOOOOOOOO 

UOrt^CiCt^O^^ON^OiSroiOOOOOOOOOOOOOOO 

o  «*pfMMNnmmm««ooooooeoooooooo 


xoiAr*oci<n<*i««t 
W«  ^  O  W) 

-  nor - - 


Kr®ro«oo«t««»T«N«in«aa'»o«A*AW«iflor©^ 


OK»«*m  I  I 
•Mill 
I 


a 

hi  *  &•  ©  •«♦  m  ^  m  »•«■  «o  o  *■• 

UU  •••••••••• 

QH  — «« 

hi  x 

u 

X 

M 

m3 

m  on«o"i9>®*nn« 

XX  OOf*®M*>®rt»MPl 

acx  o  c«  a»  «-»  *»  «■*  a*  «  — 

HH  Omm«»W 

m3  ©  ©  m  «o  «  *o  f“  r»  x  o* 

ox  o<Mm««oN®,®r 

x  . .  •  #  •  • 

X  O  M*  «4  «4  *4  «4 

o 

X 

H 

a  x 

X  hi 

•3** 

lb 

O  X 

hi 

x  m3 

hi  *3  **•*•*•*»«•*•* 

o 

m 

N 

©•»«««•<•<©•* 
ow«h««0ri<M« 
x  o  r-  ir>  m  ©  r-  «p  ^  *  r* 

x  ooo>«n«A«»® 

H  on«»ci©o«f© 

hi  . .  •  .  •  • 

o  »**<n«wrt 

X  ^fi«f©«r©©o 
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PARAMETERS  PROM  OLD  SOLUTION  ESD  RESTART  DATA 


T 


9  c*  •*  o 

9  *4 

9  «4 

•  •  1 

MMOmR 

f*  «n 

O  9 

9  9  CM  9 

9  © 

CM  •« 

oc**i# 

O  O 

CM  f" 

•  «  •  « 

«  • 

•  • 

O  CM  o 

O  O  CM 

o 

•4  9 

CM  *4 

N  *4  t*- 

9  c-  r- 

9 

1 

r-  m  *4 

9  9  9 

IM 

9  9 

CM  9  O 

CM 

u 

o  o  o 

o  o  o 

r- 

««** 

CM  O 

-M  «D 

H 

•  •  • 

•  •  • 

• 

N«OM 

9  r- 

m  *4 

X 

•  1  1 

• 

®  1 

o  r* 

9  9 

X 

O  9  O  O 

CM  *4 

O  r* 

H 

O  «4 

CM  9 

O  M 

o 

H  •* 

•  •  •  • 

•  • 

4  • 

X 

9  O  CM 

9  0  9 

o 

X 

a 

•4  CM 

CM  *4 

o 

•*  9  9 

o  r*  O 

r* 

t 

o  9  r* 

nrr 

9 

X 

9  0  9 

CM  9  O 

CM 

o 

o  o  o 

O  O  O 

r* 

0  0  9  9 

m  in 

CM  9  O 

M 

•  •  • 

•  ■  • 

• 

O 

O  O  CM  O 

r*  9 

0  9  0 

H 

•  «  1 

3  •• 

9  9  r-  9 

9  CM 

0  0  9 

X 

*  «N 

9  0  9  9 

9  9 

0  9  9 

t* 

J 

O  9  9  r> 

O  9 

O  O  f 

o 

•  •  •  • 

•  • 

•  ■  • 

9  •*  9  |- 

r»  o  9  o 

o 

«■*  H  W 

m 

9  O 

J 

0  9  0  0 

•4  CM  9  9 

9 

1 

t  9 

o 

••099 

o  o  o  o 

O 

X 

O  CM  o  9  *4 

to  «4  9  *4  O 

9 

o  o 

H 

O  O  O  O  O 

o  o  o  o  o 

o 

X  w 

M  -<  »n  n 

«  r- 

9  CM  O 

X 

X  •  •  •  • 

X  •  •  •  • 

• 

•a 

r*»  9  cm  9 

O  9 

9  9  9 

o 

o 

O  •  1  •  • 

CM 

9  o  **9 

9  9 

9  r-  o 

u 

«  •M-RM 

9  9 

•  CM  CM 

Id 

hi 

r-  «  9  o 

9  CM 

CM  CM  CM 

o 

u 

u 

•  •  •  • 

•  • 

•  •  • 

id 

X  O  CM  CM  *4 

O 

O  9 

CM  -« 9 

m 

9  r* 

o 

o 

X  •* 
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APPENDIX  C 


ATWSEQ 

CON  IN 

CPDESD 

DBNDY 

ESDCON 

ESDDRV 

FARBDY 

GCALC 

1C 

INPUT 

INTERP 

MESH 

MESH  IN 


ESD  SUBROUTINE  DESCRIPTION 


Main  program  for  both  the  ESD  and  the  FPE  code. 

Input  routine  for  CONMIN  data 

Computes  the  objective  function  and  constraints  when  optimi¬ 
zation  is  selected. 

Computes  value  of  potential  on  the  downstream  boundary 
Cycles  CONMIN  with  the  ESD  code  when  optimization  is  selected. 
Sequencing  routine  for  the  ESD  code 

Interpolates  potential  on  exterior-interior  mesh  boundaries 

Transfers  updated  values  of  circulation  to  two-dimensional 
array 

Calculates  coefficients  for  extrapolation  formulas  used  to 
find  the  jump  in  potential  at  the  wing  trailing  edge 

Main  input  routine  for  ESD  data 

interpolates  previously  stored  solution  for  use  as  initial 
guess 

Generates  stretched  x,  y  or  z  grids 

Performs  interpolation  searches  and  stores  mesh  index  locations 


MOOWNG 


OUTP 

PC  I  NT 

PFINT 

SECTIN 

SETUP 

SETUPX 

SIMP 

SLOPY 

SMTH 

SOLVE 

SOLVEX 

SPLN1 

STORE 

TCOEF 


Perturbs  wing  geometry  when  optimization  is  selected. 

Printed  output  routine 

Interpolates  potentials  onto  coarse  mesh 

Interpolates  potentials  onto  fine  mesh 

Input  routine  for  airfoil  section  coordinates 

Computes  planform  parameters  and  computational  grid  for  fine 
inner  mesh 

Computes  rectangular  coarse  exterior  mesh  and  planform 
geometry 

Simpson's  rule  numerical  integration  routine  used  to  compute 
flow  parameters 

Computes  wing  surface  slopes 

Smooths  input  section  ordinates 

Performs  SLR  solution  for  inner  grid 

Performs  SLR  solution  for  outer  grid 

Continuous  derivative  interpolation  routines 

Stores  value  of  potential  jump  in  a  single  dimension  array 

Computes  finite  difference  equation  coefficients  on  interior 
mesh 
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TCOEFX 


Computes  finite  difference  coefficients  for  coarse  exteri 
mesh 


WINGCO 

WNGBDY 


Generates  wing  boundary  conditions 

Computes  potential  at  wing  surface  by  interpolation  from 
fine  mesh  solution. 
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ATWSEQ 

8LKDAT 

CON!  N 

COORO 

CPDES 

CPUVW 

CVOL 

DINT 

EST1M 

F 

FLOW 

FORCE 

FPECON 

FPEDRV 


APPENDIX  0 

FPE  SUBROUTINE  DESCRIPTION 

Main  program  for  both  ESO  and  FPE 
Block  data  routine 
CONMIN  input  routine 

Sets  up  stretched  parabolic  and  spanwise  coordinates 

Computes  objective  function  and  constraints  when  optimization 
is  selected 

Plots  Cp  at  equal  intervals  in  the  mapped  plane 

Integrates  lift,  induced  drag,  pitching  moment  and  shock 
drag  by  control  volume  integration 

Combines  integrand  with  elemental  area  for  control  volume 
integration 

Computes  initial  estimate  of  reduced  potential 

Function  that  defines  integrand  for  control  volume  integration 

Calculates  x,  y,  z  coordinates,  u,  v,  w  velocity  components, 
Cp,  and  density  ratio  at  grid  points 

Calculates  section  force  coefficients 

Sequences  CONMIN  with  the  FPE  code 

Driver  routine  for  FPE  code 

H9 


GEOH 


Defines  geometry  of  the  wing 


GRAPH 

INPUT 

INTPL 

MDWFPE 

MIXFLO 

PPXY 

PRTMSH 

PRTPLT 

PXPYP 

RDGSO 

REFIN 

REFSMO 

RESTRT 

SAVE 

SECTIN 

SINGL 


CALCOHP  plot  routine  for  section  pressure  distributions 
Inputs  FPE  data 

Taylor  series  interpolation  routine 

Perturbs  wing  geometry  when  optimization  is  selected 

Solves  equations  for  mixed  subsonic  and  supersonic  flow 

Prints  and  plots  mesh  x,  y  values 

Prints  description  of  computational  mesh 

Sequencing  routine  for  printed  and  plotted  output 

Prints  values  of  XP  and  YP 

Transfers  current  values  of  potential  to  temporary  array 
Performs  mesh  halving 
Smooths  refined  mesh 

Reads  restart  data  for  start  from  a  stored  solution 
Saves  data  on  mass  storage  for  future  restart 
Inputs  control  station  airfoil  section  coordinates 
Generates  singular  line  for  square  root  transformation 


SHOO 


Smooths  potential  values 


SOLVE 


Sequencing  routine  for  SCR  solution 


SPUF  Spline  fit  routine 

STR2BL  Outputs  data  for  TWODBL  program 

SURF  Interpolates  mapped  wing  surface  at  mesh  points 

THREEO  Generates  three-dimensional  plots 

TOTFOR  Calculates  total  force  coefficients 

TRI  Computes  areas  and  normal  vectors  for  triangular  elements 

used  in  control  volume  integration 

VELO  Computes  surface  velocity 

WRTSGO  Transfers  potential  values  from  temporary  storage 

XYSING  Computes  coordinates  of  singular  point  at  each  span  station 

YSWEEP  Performs  row  relaxation 
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CURVFT 

FUNCT 

GRADNT 

INPUT 

INTI 

I NT2 

LAMNAR 

LGRNGE 

PRECAL 

PROF I L 


TWODBL  SUBROUTINE  DESCRIPTION 

Evaluates  a  polynomial  f(x,y)  for  a  given  set  of  coefficients 

FUNCT  expresses  the  functional  relationship  between  displacement 
thickness  6*  and  correlation  number 

Computes  the  gradient  of  a  function  with  respect  to  x 
Reads  and  prints  all  input  data 

Computes  integrand  used  in  first  call  of  SIHPS1  integration 
routine 

Computes  integrand  used  in  second  call  of  SIMPS1  integration 
routine 

Solves  the  laminar  boundary-layer  equation,  computes  laminar¬ 
boundary-layer  parameters,  checks  for  instability  and  transition 
to  turbulent  flow  and  computes  initial  values  for  turbulent 
boundary  layer  when  transition  occurs. 

Performs  four  point  Lagrange  interpolation 

Computes  initial  parameters  required  for  solution  of  boundary 
layer  equations 

Prints  all  the  principal  boundary- layer  parameters  computed  by 
LAMNAR  and  TURBLN,  and  calculates  and  prints  the  laminar  and 
turbulent  velocity  profiles 


* 
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RQOT 

RUNKUT 

SIMPS] 

SMOOTH 

SPLINE 

TURBLN 


Locates  a  root  for  a  given  function  in  a  given  interval. 

Solves  the  coupled  ordinary  differentia)  equations  of  the  turbulent 
boundary  layer  using  fourth-order  Runga-Kutta  method. 

Performs  numerical  integration  by  a  modified  Simpson's  rule 
method 

SMOOTH  is  a  simple-data  smoothing  array 
Cubic  spline  interpolation  routine 

Solves  the  turbulent-boundary  layer  equations  and  computes  the 
turbulent-boundary  layer  parameters 
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WING-PYLON/NACELLE  PROGRAM  SUBROUTINE  DESCRIPTION* 

AXIS  Computes  finite  difference  terms  along  nacelle  axis 

CUTOUT  Aborts  program  when  error  condition  is  detected 

GRID  Generates  finite  difference  mesh  for  nacelle 

GRIDGEU  Generates  a  wing  grid  tailored  to  the  pylon/nacelle  combination 

MAIN  Computes  transonic  flow  about  wing 

METRIC  Computes  metric  coefficients  for  finite  difference  equations  in 

nacel le  solution 

NACELLE  Computes  transonic  flow  about  nacelle 

OUTNAC  Outputs  nacelle  potential  solution  to  Unit  9 

PYLCO  Computes  pylon  first  and  second  derivatives 

PYLON  Computes  special  diagonal  and  RHS  difference  expressions  for 

wing-pylon  solution. 

SETIN  Computes  coefficients  of  the  difference  equations  for  tridiagonal 

solver  in  nacelle  solution. 

SHAZ  Transfers  wing  and  nacelle  boundary  conditions  during  wing/nacelle 

solution  cycle 

*  See  Appendix  C  for  description  of  subroutines  used  by  ESD  wing  code  but  not 
described  here. 
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SLOPPY 

TALA 

TRIO 


Sets  up  computation  of  wing  surface  slopes 


Main  program  for  wing-pylon/nacelle  code 


Tridiagonal  equation  solver 
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TAPE4 

TAPE5 

TAPE6 

TAPE7 


TAPE8 


TAPE9 

TAPEJJO 


LOGICAL  UNITS  REFERENCED  IN  PROGRAMS 

Input  unit  for  stored  solution  data  used  to  restart  FPE  program 
Card  reader  logical  unit  (all  programs) 

Printer  logical  unit  (all  programs) 

Output  unit  for  solution  data  used  to  restart  subsequent  runs  of 
the  FPE  program 

Input  unit  in  TW006L  program  for  data  generated  in  design  and 
analysis  programs 

Input  unit  in  LINKUP  and  LINKDN  for  data  generated  by  TWODBL 

Input  unit  for  airfoil  section  data  when  AFIN-F  in  ESD  and  FPE 
programs 

Input  unit  for  solution  data  used  to  restart  nacelle  solution 
in  wing-pylon/nacel le  code 

Output  unit  for  wing  ordinates  generated  by  LINKUP  and  LINKDN 

Output  unit  for  TVODBL  data  used  in  LINKUP  and  LINKDN 

Output  unit  for  nacelle  solution  data  used  to  restart  subsequent 
runs  of  the  wing-pylon/nacelle  code 

Input  unit  for  stored  wing  solution  data  used  to  restart  the 
ESD  and  wing-pylon/nacelle  codes 
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TAPE]! 


TA PE]4 


Output  unit  for  wing  solution  data  to  be  used  to  restart  subse¬ 
quent  runs  of  the  ESD  and  wing-pY ton/nacel le  codes. 

Output  unit  in  ESD  and  FPE  codes  for  data  to  be  used  in  TWODBL 
program 

Input  unit  in  wing-pylon/nacelle  code  for  airfoil  section  data 
when  AFIN«F. 

ESO  plot  data  output  unit 
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